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ABSTRACT 


Study  is  made  of  the  vibrational  behavior  of  flat  plates  when 
placed  parallel  to  a  uniform  stream  and  allowed  one  degree  of  free¬ 
dom  torsional  oscillation  around  a  vertical  axis  through  their  leading 
edge. 


The  experimental  work  was  carried  out  in  the  7-1/2" x  9 "  test 
section  of  an  open-circuit  water  tunnel  in  the  Hydrodynamics  Laboratory 
of  the  Civil  Engineering  Department  at  M.I.T.  Four  test  plates  were 
used,  each  1/h  inch  thick  and  2  inches  in  chord,  but  having  different 
trailing  edge  geometries.  Measurements  were  made  of  the  amplitude  and 
frequency  of  plate  vibration  and  of  the  fluctuating  total  head  in  the 
wake  over  the  range  of  free-stream  velocity  producing  significant 
vibrational  excitation.  Measurable  vibrations  are  found  to  occur  in 
three  distinct  regions  surrounding  the  natural  resonant  frequencies  of 
the  system. 

An  equation  of  motion  of  the  plate-spring  system  is  developed 
incorporating  the  hydrodynamic  loads  given  by  the  linearized  potential 
theory,  and  the  unknown,  vortex-induced,  forcing  moments .  Using  this 
linear  model  and  the  experimental  data,  forcing  moment  coefficients  are 
obtained  which  show  a  unique  (for  all  plates)  and  piecewise  linear  de¬ 
pendence  upon  a  "self -excitation"  parameter  consisting  of  the  ratio  of 
the  mean  transverse  velocity  of  the  trailing  edge  to  the  free-stream 
velocity.  This  dependence  of  moment  coefficient  on  vibrational  ampli¬ 
tude  is  equivalent  to  a  negative  damping  and  thus  yields  eigenvalues 
of  frequency  at  which  the  amplitude  tends  to  grow  without  bound. 

These  critical  frequencies  and  the  free-stream  velocities  at 
which  they  occur  are  predicted  knowing  only  the  physical  system  con¬ 
stants,  the  body’s  steady-state  Strouhal  number,  and  the  slope  of  the 
moment  coefficient  function.  They  are  shown  to  agree  with  the  bounds 
of  the  region  of  large  amplitude  motion  and  variable  Strouhal  number 
obtained  in  these  experiments  and  in  those  of  other  investigators. 

Major  conclusions  of  this  study  are: 

1.  Controlled  variation  of  the  system  parameters,  i.e.,  in¬ 
creasing  the  combined  structural-viscous  damping  and  increasing  the 
steady-state  Strouhal  number,  allows  design  of  a  system  in  which  self¬ 
excitation  is  avoided. 

2.  The  primary  role  of  trailing  edge  geometry  in  flow-induced 
vibrations  is  in  determining  the  free-stream  velocity  for  resonant 
response  through  establishment  of  the  steady-state  Strouhal  number. 

3.  Self -excitation,  the  increase  in  moment  coefficient  with 
increasing  vibrational  amplitude,  is  due  both  to  increasing  two- 
dimensionality  of  the  wake  and  to  increasing  strength  of  the  discrete 
vortices . 
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I .  INTRODUCTION 


A.  Flow-Induced  Vibrations 


It  is  an  observed  fact  that  structural  members  subjected  to  a 
flow  field  may,  under  certain  conditions,  suffer  considerable  vibrations. 
Vibration  of  a  structure  is  undesirable  for  many  reasons,  one  of  which 
is  the  danger  of  collapse  from  fatigue  of  the  material.  It  is  the  aim 
of  the  structural  designer  to  avoid  such  a  possibility. 

Vibrations  of  smokestacks,  pipe-lines,  telephone  wires,  submarine 
periscopes,  television  antennas,  turbine  and  propeller  blades,  vanes  and 
the  flutter  of  aircraft  wings,  are  but  a  few  examples  of  flow- induced 
vibrations . 

This  phenomenon  has  long  attracted  the  attention  of  engineers 
and  scientists  and  is  still  not  completely  understood,  remaining  a 
stimulating  field  of  analytical  and  experimental  research.  From  the 
fluid  mechanics  point  of  view  structures  can  be  classified  into  two  cate¬ 
gories:  1.  Streamlined  structures  in  which  the  fluid  follows  continu¬ 
ously  the  geometric  contour  of  the  body  and  2.  Unstreamlined  or  "bluff" 
structures  in  which  the  flow  becomes  separated  from  the  boundary. 

Structures  of  the  first  category  are  used  extensively  in  aerodynamics 
and  well  established  theories  exist  for  the  computation  of  their  dynamic 
loads.  Structures  of  the  second  category  are  used  extensively  in  civil 
and  mechanical  engineering  practice.  Separation  of  flow  occurs  in  such 
structures  and  since  there  is  no  complete  theory  for  computation  of  the 
resulting  dynamic  loads,  experimental  work  becomes  the  basic  source  of 
information. 

The  shedding  of  discrete  vortices  in  the  wake  of  an  unstream¬ 
lined  body  submerged  in  a  uniform  stream  has  been  recognized  as  instru¬ 
mental  in  creating  transverse  periodic  forces  on  the  body.  The  periodic 
discharge  of  discrete  vortices  in  the  wake  gives  rise  to  periodic  circu¬ 
lation  around  the  body  and  thus  to  a  periodic  lift  force,  the  magnitude 
of  which  is  proportional  to  the  strength  of  the  discrete  vortices.  As 
the  frequency  of  periodic  lift  forces  approaches  a  natural  vibrational 
mode  of  the  structure  substantial  vibrations  may  occur.  It  is  thus  obvious 
that  a  detailed  knowledge  of  the  flow  field  in  the  wake  of  the  body  is  of 
paramount  importance  in  estimating  the  frequency  and  magnitude  of  the  perio¬ 
dic  lift  forces. 

B.  Scope  of  the  Investigation 

The  investigation  reported  herein  is  part  of  a  fundamental  re¬ 
search  program  on  wake  mechanics  undertaken  in  the  Hydrodynamics  Labora¬ 
tory,  Department  of  Civil  Engineering,  at  MIT  under  the  sponsorship  of 
the  David  Taylor  Model  Basin,  Office  of  Naval  Research.  It  is  the  general 
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aim  of  this  project  to  contribute  to  the  understanding  of  the  problem  of 
flow-induced  vibrations.  The  flat  plate  has  been  selected  as  the  struc¬ 
tural  element  to  be  tested  since  it  constitutes  the  most  simple  structural 
element  similar  in  form  to  the  blades  and  vanes  of  primary  Interest  to  the 
Naval  hydrodynamicist.  Previous  investigators  have  found  the  trailing 
edge  geometry  to  be  of  tremendous  importance  in  the  vibrational  problem, 
which  accounts  for  its  selection  as  a  primary  variable  in  the  present  tests. 

In  order  to  simplify  the  problem  and  to  reduce  the  important 
factors  to  a  minimum,  the  plates  are  tested  at  zero  mean  angle  of  attack 
(parallel  to  the  ambient  flow) .  The  plates  are  provided  with  an  external 
elastic  restoring  force  and  the  motion  of  the  plate  is  restricted  to  one 
degree  of  freedom,  that  of  oscillatory  motion  around  an  axis  containing 
the  leading  edge. 

The  first  phase  of  this  project  dealt  with  an  exploratory  experi¬ 
mental  investigation  of  the  vibrational  behavior  of  the  plates  as  related 
to  trailing  edge  geometry,  externally  supplied  restraint  and  mean  ambient 
flow  velocity.  A  theoretical  review  of  the  vibrational  problem  and  com¬ 
parison  of  the  experimental  results  with  the  theoretical  concepts  closed 
this  phase  (1) . 

The  second  phase  dealt  with  an  experimental  investigation  of  the 
early  wake  structure  of  the  plates  when  vibrational  motion  was  prevented 

(2)(3)<U). 

It  is  the  aim  of  the  present  investigation  to  provide  a  quantita¬ 
tive  estimate  of  the  periodic  forces  acting  on  the  plate  and  to  correlate 
them  with  pertinent  wake  characteristics.  For  this  purpose  provision  is 
made  for  experimental  determination  of  the  structural  characteristics 
(natural  frequency  and  damping)  of  the  plate-spring  system.  Amplitudes  and 
frequencies  of  vibration  are  measured  to  establish  experimentally  the 
response  curve  of  every  plate-spring  system  throughout  the  velocity  range 
of  measurable  vibrational  excitation.  Root-mean-square  and  spectral 
measurements  of  the  power  density  of  the  total  head  fluctuations  are  made 
in  the  early  wake. 


II.  LITERATURE  REVIEW 

Fairly  exhaustive  surveys  of  previous  work  in  the  areas  of  flow- 
induced  vibrations  and  wake  mechanics  are  given  elsewhere  (1) (2) (U) ,  thus 
only  a  brief  summary  is  included  here. 

A.  Engineering  Aspects  of  Flow-Induced  Vibrations 

The  vast  majority  of  studies  have  treated  either  the  circular 
cylinder  or  the  thin  flat  plate  and  have  centered  around  one  or  more  of 
the  following  questions: 
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1*  What  Is  the  origin  of  the  exciting  force? 

2.  Are  the  large  amplitude  vibrations  self-excited  or  do  they  represent 
a  forced  vibration  at  a  resonant  structural  frequency? 

3.  What  means  can  be  employed  to  prevent  large  amplitude  vibrations? 

Dockstader  (5)  attributes  the  vibrations  of  steel  stacks  caused  by  wind  to 
forces  induced  by  the  von  Karman  vortices.  Parmakian  and  Jacobson  (6)  re¬ 
porting  on  vibrations  of  a  slow-speed  hydraulic  turbine  attribute  the 
vibrations  to  vortex  shedding  at  the  trailing  edge  of  the  blades. 

Donaldsen  (7),  Heskeetadand  Olberts  (8),  Ippen,  et  al.  (1),  point 
out  the  importance  of  the  trailing  edge  geometry  on  the  maximum  amplitudes 
obtained  in  tests  with  flat  plates  and  consider  the  vibrations  as  induced 
by  the  vortex  shedding. 

The  singing  of  propellers  is  attributed  by  Conn  (9)  to  torsional 
vibrations  of  the  blades  caused  by  hydrodynamic  forces  which  change  in 
magnitude,  direction,  and  location.  Hughes  (10)  attributes  the  singing 
of  propellers  to  a  rapid  change  of  forces  as  the  propeller  blade  passes 
through  the  variable  wake.  Gongwer  (11),  from  experimental  studies  on 
vanes  singing  in  water,  reaches  the  conclusion  that  singing  is  caused  by 
the  action  of  the  von  Karman  vortex  street. 

Dickey  and  Woodruff  (12)  reporting  on  field  data  of  steel  stack 
vibrations  consider  the  important  question  to  be  whether  the  vibrations 
are  of  the  self -excited  type  (negative  air  damping)  or  forced  vibrations 
at  a  structural  resonance.  Strouhal  numbers  computed  from  the  field  data 
differ  appreciably  from  the  experimentally  determined  values  for  rigid 
cylinders,  being  in  general  smaller. 

Penzien  (13)  from  experiments  in  a  wind  tunnel  with  cylinders  con¬ 
cludes  that  the  vibration  of  tall  stacks  are  self -excited  rather  than 
forced.  Marris  (lit)  from  experiments  on  canitlever-supported  probes  con¬ 
cludes  that  negative  hydrodynamic  damping  is  responsible  for  the  large 
amplitude  vibrations  in  the  vicinity  of  structural  resonance. 

From  the  above  brief  survey  of  the  literature  it  can  be  deduced 
that  while  the  shedding  of  vortices  in  the  wake  is  recognized  as  the 
initiating  mechanism  of  what  is  truly  a  forced  vibration,  other  effects, 
such  as  negative  hydrodynamic  or  aerodynamic  damping,  may  become  important 
as  the  amplitude  of  vibration  grows.  Under  these  latter  conditions  the 
motion  may  be  considered  self-excited  with  the  amplitude  growing  continu¬ 
ously  until  limited  by  an  increase  in  the  damping  or  until  structural 
failure  occurs. 

Means  of  preventing  flow- induced  vibrations  have  been  studied. 

The  enclosure  of  a  cylinder  in  a  concentric  perforated  shroud  has  been 
found  as  an  effective  suppressor  of  vibrations  by  Price  (15)  after  exten¬ 
sive  tests  in  a  wind  tunnel.  Design  criteria  for  helical  aerodynamic 
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spoilers  have  been  formulated  by  Weaver  (16)  for  cylinders .  Baird  (17) 
reporting  on  the  vibration  of  a  pipe-line  suspension  bridge  excited  by 
wind  indicates  that  the  use  of  a  sawtooth  paneling  around  the  pipe-line 
proved  very  satisfactory  in  minimising  the  vibrations.  Reduction  of  the 
thickness  of  the  trailing  edge  of  the  blades  proved  satisfactory  in  re¬ 
ducing  the  vibrations  of  turbines  (6) .  Sharp  and  thin  leading  edge  of 
the  propeller  blades  is  suggested  as  a  remedy  for  the  singing  of  propel¬ 
lers  by  Shannon  (18)  and  Hughes  (10) . 

B.  Wake  Studies 


The  shedding  of  vortices  in  the  wake  of  a  body  submerged  in  a  uni¬ 
form  stream  was  understood  to  be  the  initiating  mechanism  of  vibrations 
in  many  cases  of  flow- induced  vibrations.  Thus  the  wake  studies  and  the 
detailed  description  of  the  flow  field  in  the  wake  is  of  basic  importance 
in  estimating  the  vortex  induced  forces  on  a  structure. 

A  thorough  survey  of  wake  flows  for  cylinders  and  flat  plates  is 
given  in  References  (1) (2) (ii) .  In  the  present  brief  summary  the  aspects 
of  wake  flow  related  to  the  vibrational  problem  will  be  stressed.  It 
should  be  stated  that  a  complete  theoretical  model  for  the  wake  flows 
does  not  exist  and  we  must  rely  heavily  upon  experimental  information. 

The  first  mathematical  model  for  the  description  of  the  wake  of  a 
bluff  body  was  given  by  Kirchhoff  (19)  in  1869  and  was  a  big  advance  in 
classical  hydrodynamics,  removing  the  obstacle  of  D’Alemberts  paradox. 

Based  on  this  model  and  using  the  momentum  principle,  the  drag  coefficient 
on  the  body  can  be  obtained.  Many  of  the  basic  assumptions  of  this  model 
are  not  substantiated  experimentally  and  drag  computations  underestimate 
the  true  experimental  values. 

The  second  mathematical  model  is  the  one  presented  by  von  Karman 
(20) (21)  in  19H-12.  It  represents  the  wake  by  two  parallel  rows  of 
equally  spaced  and  staggered  point  vortices  in  an  otherwise  irrotational 
field  of  flow}  the  important  parameters  being  the  longitudinal  and  trans¬ 
verse  spacing  of  vortices  and  the  strength  of  the  individual  vortex.  Von 
Karman  proved  that  the  only  stable  configuration  of  this  form  is  the  one 
with  a  spacing  ratio  h/a  ■  0.281.  Rosenhead  (22)  substantiated  the  above 
stability  ratio  and  proved  that  the  configuration  possesses  instabilities 
of  higher  order.  Most  of  the  experimental  work  done  on  wake  flows  is 
centered  on  the  above  model  and  provides  data  both  supporting  and  rejecting 
and  applicability  of  the  model. 

The  von  Karman  model  provides  adequate  explanation  of  the  periodic 
nature  of  the  wake  which  is  instrumental  in  creating  periodic  transverse 
forces  on  the  body. 

Theoretical  and  experimental  efforts  were  made  in  two  directions; 
in  establishing  the  actual  spacing  ratio  h/a  of  the  vortex  street  and  in 
predicting  the  strength  of  the  individual  vortex. 


Experimental  work  (23)  (2ii)  (25)  (26)  (2?)  has  proved  that  the  trans¬ 
verse  spacing  of  the  vortex  street  increases  downstream  in  the  wake  and 
Birkhoff  (28)  justifies  this  widening  on  mathematical  grounds.  Efforts 
were  made  to  connect  the  transverse  spacing  of  the  vortices  to  the  body- 
thickness. 

Heisenberg  (29)  and  Prandtl  (30)  attempted  to  connect  the  strength 
of  the  individual  vortex  with  the  rate  of  vorticity  discharge  into  the 
wake  from  each  side  of  the  body  by  the  separated  shearing  layers. 

The  Strouhal  number,  a  parameter  of  importance  in  the  description 
of  periodic  wakes,  has  received  special  attention.  For  a  certain  body 
geometry  it  has  been  established  that  the  Strouhal  number  is  a  function 
of  body  Reynolds  number.  Reported  experimental  values  of  Strouhal  number 
are  given  for  different  body  geometries  and  Reynolds  numbers  in  References 
(31) (32) (33) (3U) ,  (1) (3) (3),(ll) .  The  variety  of  these  experimental  values 
reveals  the  importance  of  body  geometry  and  their  range  for  the  same  body 
studied  by  different  investigators  suggests  a  sensitivity  to  experimental 
method. 


The  wake  flow  is  usually  classified  by  the  body  Reynolds  number. 

A  thorough  experimental  work  for  cylinders  is  reported  by  Roshko  (31)  and 
sufficient  information  for  the  flat  plate  by  Taneda  (25>) ,  While  the 
shedding  of  vortices  gives  an  easy  explanation  of  periodic  transverse 
forces  acting  on  the  body,  rigorous  theoretical  prediction  of  these  forces 
has  not  yet  been  developed.  Actual  direct  measurements  of  vortex  induced 
lift  forces  have  been  reported  (35) (36) (37)  only  for  cylinders  in  spite  of 
the  importance  in  design  of  other,  more  streamlined,  body  shapes. 

Most  of  the  theoretical  and  experimental  efforts  toward  a  des¬ 
cription  of  the  wake  have  been  centered  around  two-dimensional  models 
while  in  actuality,  the  total  vortex  induced  lift  forces  acting  on  a  body 
are  integral  quantities  over  the  span  of  the  body.  Thus  three-dimensional 
effects  can  greatly  influence  their  value.  Thorough  investigations  of  the 
three-dimensional  structure  of  wakes  do  not  exist.  Roshko  (31)  reports 
some  experimental  observations  on  cylinders. 

Thus  far  the  above  discussion  has  been  based  on  wakes  of 
stationary  bodies,  however,  the  lift  forces  associated  with  the  shedding 
of  vortices  can  cause  body  vibrations  which  may  greatly  influence  the 
wake  structure. 

The  problem  of  the  thin  flat  plate  oscillating  in  a  uniform  in¬ 
compressible  flow  has  been  solved  by  Theodorsen  (38)  on  the  basis  of 
linearized  potential  theory  and  this  solution  constitutes  the  basis  of 
flutter  analysis  in  aerodynamics.  For  bluff  bodies  where  separation  and 
shedding  of  vortices  occurs  such  a  solution  does  not  exist. 
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A.  Water  Trawl 


III.  EXPERIMENTAL  EQUIPMENT 


The  experiments  reported  herein  were  performed  in  a  closed- jet, 
open  circuit  water  tunnel  which  has  been  described  in  detail  in 
References  (l) (39) •  A  schematic  view  is  presented  in  Fig.  1. 

An  effort  was  made  to  minimize  the  free  turbulence  level  in  the 
test  section.  It  was  proved  that  by  using  pump  No.  Ill  (rated  capacity 
3000  GPM  at  H  ■  and  1175  RPM)  and  running  the  water  through  the 
longest  approach,  path,  a  minimum  free  turbulence  level  in  the  test  section 

is  obtained  (\/u'2  /V  -  3?o  or  h?0) . 

The  uniformity  of  the  flow  in  the  test  section  was  checked  and 
the  tunnel  contraction  was  recalibrated  for  flow  measuring  purposes. 

B.  Test  Plates 


Four  flat  plates  with  different  trailing  edge  geometries  were 
used.  The  cross  section  of  the  plates  and  the  trailing  edge  configura¬ 
tions  are  shown  in  Fig.  2.  The  plates  are  mounted  vertically  in  the  test 
section,  parallel  to  the  flow  in  such  a  way  as  to  permit  rotation  about 
the  center  of  their  semi-circular  leading  edge.  They  are  fabricated  of 
aluminum  and  all  but  plate  D  have  been  electroplated  to  prevent  corrosion. 
Plate  D  has  a  somewhat  rougher  surface  than  the  three  others  due  to 
corrosion  resulting  from  its  use  in  a  previous  investigation  (2) . 

C.  Torsion  Spring  and  Plate  Mounting  Device 

1.  Basic  Design  Criteria 

The  experimental  results  reported  in  Reference  (1)  served  as 
a  guide  in  the  design  of  the  spring  and  plate  mounting  assembly.  The 
criteria  adopted  were  the  following: 

a)  Single  degree  of  freedom  motion  of  the  plate.  Oscilla¬ 
tory  motion  of  the  plate  as  a  rigid  body  around  an  axis  through  the 
center  of  the  semi-circular  leading  edge  will  be  permitted.  No  other 
motion  should  be  permitted  in  order  to  avoid  coupled  oscillations. 

b)  The  spring  should  be  sensitive  enough  to  provide  for 
accurate  measurement  of  very  small  angular  displacements  of  the  plate. 

c)  The  plate-spring  system  and  the  whole  mounting  assembly 
should  provide  for  the  minimum  possible  structural  damping. 

d)  Measurement  of  moments  should  be  made  on  the  flow  side  of 
any  water  seals. 
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CIRCULATION  PUMPS  and  DISTRIBUTION  CIRCUIT 

( PLAN  ) 


IJ '  n  **4  in 

Pump  I 

Fig.  1  Water  Tunnel  and  Distribution  Circuit 
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PLATE  A  - ► 

PLATE  B  - - 

PLATE  C  *- 

PLATE  D  — 


Fig.  2  Test  Plates 


2.  Description  of  the  Mounting  Assembly 

The  plate  mounting  device  is  shown  in  Figs.  3,  U>  and 

Two  rectangular  slots  on  the  upper  and  lower  lucite  faces  of  the 
test  section  permit  the  plate  to  be  inserted  vertically.  A  cylindrical 
shaft  is  clamped  tightly  on  each  end  of  the  plate,  the  axis  of  the  shafts 
coinciding  with  the  axis  of  rotation  previously  specified.  Each  shaft 
carries  a  stainless  steel,  doubly  shielded,  radial  ball  bearing,  which 
fits  a  counterbored  hole  in  an  aluminum  plate  attached  to  the  wall  of  the 
test  section.  Sealing  of  both  plates  is  provided  by  O-rings  and  their 
exact  positioning  is  insured  by  keying.  A  flanged  aluminum  pipe  is 
pinned  and  fastened  to  the  upper  aluminum  plate.  The  top  of  the  pipe  is 
covered  with  an  aluminum  plate  carrying  a  central  hole  for  insertion  of 
the  torsion  spring. 

The  spring  consists  of  a  hollow  cylindrical  stainless  steel 
tube  1/2"  OD  and  one  foot  long.  The  first  spring  (I)  used  had  a  wall 
thickness  of  0.035"  and  the  second  one  (II),  0.0U9".  At  its  lower  end 
the  tube  is  rigidly  fixed  to  a  cylindrical  aluminum  piece  which  carries 
a  square  hole  to  receive  the  end  of  the  upper  test  plate  shaft.  At  its 
upper  end  the  spring  tube  is  welded  to  a  stainless  steel  plate  which  is 
fastened  to  the  cover  of  the  aluminum  pipe.  The  aluminum  pipe  has  two 
small  circular  openings  to  facilitate  the  placement  of  screws  damping 
the  spring  to  the  shaft  and  the  inplace  calibration  of  the  spring. 

The  whole  mounting  device  is  assembled  in  place  after  setting 
the  plate  at  zero  angle  of  attack  (±0.001  rad) .  To  prevent  leakage 
through  the  lower  shaft-bearing  a  sealed  cover  plate  is  used.  The  whole 
upper  device  is  enclosed  in  a  lucite  casing.  Electrical  leads  from  the 
torsion  spring  strain  gages  are  passed  through  a  sealed  hole  in  this  casing. 
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All  cavities  in  the  test  section  walls  are  plugged  with  specially 
made  lucite  pieces  in  order  to  avoid  major  disturbances  of  the  flow.  The 
shoes  described  in  Reference  (It)  were  used  to  minimize  end  effects  (see 
Fig.  6). 

D.  Measurement  of  Plate  Motion 


1.  Strain  Gages 

Tatnall  type,  C9-121-R2C,  strain  gages  were  cemented  and  water¬ 
proofed  on  the  surface  of  the  torsion  springs  in  the  conventional  array 
for  torsion  measurement.  The  four  gages  form  a  Wheatstone  bridge  circuit 
supplied  by  a  12  Volt  automobile  battery,  A  variable  resistance  in  the 
circuit  provides  for  balancing  of  the  bridge.  The  voltage  output  of  the 
bridge  is  proportional  to  the  angle  of  rotation  of  the  spring. (See 
Appendix  B) . 

2.  Amplifier 

The  voltage  output  of  the  Wheatstone  bridge  is  amplified  with 
a  KIN-TEL  DC  Amplifier .  This  amplifier  has  a  flat  response  over  the 
operating  range  of  frequencies  (0  to  750  cps) .  The  design  gain  accuracy 
is  ±  l/o  from  DC  to  2KC. 

3.  Oscilloscope 

A  Dumont  dual-beam  oscilloscope  was  used  during  the  tests  for 
visual  observations  of  the  vibrational  signal. 

1*.  AC  Voltmeter 


The  total  output  of  the  Wheatstone  bridge  was  measured  with  a 
Balantine  vacuum  tube  voltmeter,  which  is  a  true  RMS  indicating  meter 
for  sinusoidal  wave  forms. 

5>.  Vibration  Analyzer 

The  analyzer  used  was  a  General  Radio  Type  7&2-B  Vibration  ana¬ 
lyzer.  It  is  a  constant  percentage  band  width  tunable  voltmeter  with  a 
logarithmic  scale.  Its  operating  range  of  frequencies  is  from  2.5  to 
750  cps.  This  instrument  is  battery  supplied,  therefore  frequent  calibra¬ 
tions  were  made.  Calibration  procedure  is  described  in  Reference  (1*3) •  A 
damping  capacitance  was  inserted  in  parallel  with  the  analyzer  meter  to 
reduce  needle  fluctuations.  Vibrational  frequencies  as  well  as  amplitudes 
of  every  component  were  measured.  The  electrical  components  for  vibrational 
measurements  are  shown  in  Figs.  1*  and  7.  A  schematic  block  diagram  of  in¬ 
strumentation  for  vibrational  measurements  is  shown  in  Fig.  8. 
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E.  Velocity  Measurements 

The  mean  velocity  in  the  test  section  (ambient  velocity)  was 
measured  using  the  water  tunnel  contraction  as  a  flow  meter  (1) .  A 
differential,  manometer  with  Meriam  fluid  No.  D-8325  (specific  gravity 
1.75)  And  water  was  used.  Mean  local  velocities  in  the  wake  of  the 
plate  were  measured  with  the  total  head  and  static  pressure  tubes  des¬ 
cribed  in  Reference  (2)  (see  Fig.  9) ,  using  open  mercury  manometers. 

F.  Measurements  of  the  Total  Head  Fluctuations 


The  gage  for  measuring  the  total  head  fluctuations  has  been 
described  in  detail  in  Reference  (2)  and  (1*0) ,  along  with  the  elec¬ 
trical  components  required.  The  measuring  system  consists  of: 


1.  instantaneous  total  head  gage, 

2.  electrometer  preamplifier, 

3.  Kin- Tel  amplifier, 

1*.  oscilloscope  for  visual  observation  of  the  signal, 

£.  vibration  analyzer  for  spectral  analysis, 

6.  thermocouple  voltmeter  which  is  a  true  RMS  indicating  meter 
for  a  signal  of  ary  wave  form. 


The  electrical  components  for  total  head  fluctuation  measurements  are 
shown  in  Fig.  10.  A  schematic  block  diagram  of  the  above  components  is 
shown  in  Fig.  11.  Calibration  of  the  gage  is  reported  in  References  (2) 
and  (3) .  The  only  difference  in  the  calibration  procedure  followed  in 
the  present  tests  was  that  the  Gage-Preamplifier- Amplifier  (at  Ga  -  £0) 
were  connected  together.  The  gage  was  calibrated  at  the  beginning  and  at 
the  end  of  the  tests. 


The  gage  sensitivity  obtained  in  the  present  tests  and  used  in 
the  analysis  of  the  data  taken  is: 


KG 


lit.60 


mv 

ft .water 


The  preamplifier  gain  was  also  checked  and  found  to  be 


■0.350  g 

mv 
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Seal*  (Inch**)  ELEVATION  VIEW 

Fig.  9  Traversing  System 
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Upstream  Edge  of  Test  Section 


17.  FORMULATION  OF  THE  LINEARIZED  VIBRATIONAL  PROBLEM 


The  plate-spring  system  employed  in  the  present  tests  was  designed 
to  insure  a  single  degree  of  freedom  torsional  oscillation  about  the 
leading  edge  of  the  plate. 

Assuming  linear  damping ,  the  equation  of  motion  for  this  plate¬ 
spring  system  can  thus  be  written  in  the  form: 

V  +  V  +  Ksa  •  Mt  (1) 


in  which: 

e  •• 

a,  a,  a  represent  the  angular  plate  displacement,  velocity, 
and  acceleration  respectively  (in  radians,  rad.  per 
sec.  and  rad/sec2) 

I  *  polar  moment  of  inertia  of  the  plate  about  the 
rotational  axis, 

c  ■  structural  damping  coefficient, 
s 

K  *  structural  spring  coefficient, 
s 

and  is  the  total  hydrodynamic  moment  acting  on  the  plate. 

The  moment  may  be  thought  of  as  consisting  of  several  parts: 

^t  theodorsen  +  ^viscous  +  ^forcing  ^ 

Theodorsen  Moment 


The  Theodorsen  moment  (38)  is  that  moment  due  to  the  potential-flow 
pressure  distribution  about  the  oscillating  plate.  For  thin  flat  plates 
undergoing  small  amplitude,  torsional,  simple  harmonic  oscillations  about 
the  leading  edge  in  a  uniform  stream  of  incompressible  fluid,  Bisplinghoff 
et  al.  (i*l)  show  this  moment  to  be: 

theodorsen  **  -1.12^np  b^  H'4 

-  |  npb3  H  V  (1  +  C(I0)  4 

-  npb2  H  V2  C(K)  a  (3) 
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in  which 


b  ■  half  chord  of  plate  (1"  in  the  present  tests) 
H  -  plate  span  (8.625"  in  the  present  tests) 

V  -  ambient  free  stream  velocity 
C(X)  ■  Theodorsen  function 
C(X)  -  F(X)  +  i  G(X)  (see  Fig.  12) 

and 


(3) 

(4) 


where  f0  is  the  frequency  of  plate  oscillation  in  cps.  For  the  present 
plate  geometry  (2b  ■  2"  and  t  ■  1/4"),  b  ■  4  plate  thicknesses  or  b  ■  4t. 
Thus  utilizing  the  plate  thickness  Strouhal  number: 


V 


(5) 


we  find: 


X  ■  8nS^. 


(6) 


Using  the  subscript  "p"  for  "potential"  we  may  write  the  Theodorsen  moment: 


in  which: 


^Theodorsen  **  M  ■-ilx-ca-Ka 
P  P  P  P 


(7) 


and 


Inertia  coefficient  (real) 

I  -  1.125  np  b^H 

Damping  coefficient  (complex) 

cp  -  \  np  b3H  V  (1  +  C(X)) 

Spring  coefficient  (complex) 

Kp  -  npb2H  V2C(X) 


(8) 

(9) 

(10) 


For  the  range  of  Strouhal  numbers  obtained  in  the  present  tests  K>  2.5. 
Under  these  conditions  the  real  part  of  the  Theodorsen  function  C(X) 
tends  to  the  value  of  1/2  while  the  imaginary  part  becomes  very  small. 
We  will  thus  neglect  the  imaginary  part  whereupon,  C(X)  «  1/2. 
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1.0 


In  this  case  Eqs.  (9)  and  (10)  became: 

Cp  »  ^  np  b^  H7 

(Si) 

Kp  -  \  np  b2HV2 

(10a) 

Viscous  Moment 


The  viscous  moment  is  a  damping  moment  due  to  the  finite  fluid 
viscosity.  It  is  actually  related  to  the  relative  fluid-plate  velocity 
but  will  be  assumed  here  to  be  given  by  the  linear  relationship: 

Won.  ■  M,  •  -V  <u> 


Forcing  Moment 

This  is  the  moment  which  must  be  present  in  order  to  initiate  plate 
motion.  It  arises  from  the  fluctuations  in  pressure  distribution  on  the 
plate  resulting  from  flow  separation  at  the  trailing  edge  and  hence  will 
likely  be  modified  by  the  magnitude  of  the  motion  which  it  initiates. 


Let  Mforcjjjg  “  M.  Equation  (l)  can  now  be  written: 


(Is+Ip)«  +  (C8+Vcp)“  +  (Ks+Kp)a  "  M  (12) 

From  purely  dimensional  considerations,  the  simplest  form  of  the  forcing 
moment  is: 


M  -  mf  |  pH(2b)2V2 


(13) 


in  which  m^  *  forcing  moment  coefficient.  When  the  forcing  moment  is 
periodic  as  in  the  case  of  excitation  by  vortex  shedding,  the  moment  co¬ 
efficient  is  conveniently  written  in  the  complex  form: 


m « 


"fo  W 


(111) 


in  which: 

f  ■  frequency  of  vortex  shedding  from  one  side  of  the  plate 
v  in  cps , 

0  -  phase  angle  between  forcing  moment  and  the  motion  a, 

and 

t  “  time  in  seconds. 
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Equation  (12)  is  then  written: 

(I  +1  )V  +  (c  +c  +c  )a  +  (K  +K  )a 
'  S  p'  '  s  v  p  v  s  p 


m 


To  2 


1  pH(2b)V 


(15) 


Under  conditions  of  steady  ambient  flow  and  small  amplitude  simple  harmonic 
vibrations,  the  coefficients,  Is,  Ip,  cs,  cv,  Cp,  Kg,  and  Kp  may  be  expected 
to  be  independent  of  time  and  plate  motion.  Assuming  mfQ  to  be  independent 
of  time,  equation  (15)  represents  the  damped  forced  vibration  of  a  linear 
single  degree  of  freedom  system.  The  solution  of  the  equation  (15)  is 
well  known.  The  vibrational  frequency  is  equal  to  the  forcing  frequency, 
i  .6 . 


f  =  f 
a  v 


and 


Letting: 


i2nft 


“o  6 


I  +  I 
s  p 


c  +  c  +  c 

S  V  p 


(16) 

(17) 

(18a) 

(18b) 


K.  -  K  +  K  , 
t  s  p  ’ 


the  amplitude  response  is  given  by: 

1 

nfo  2 


m.  I  pH(2b)2V2 


r T  2,  2  2v2  2  2^/2 

[It  (con  -CO  )  +  ct  CO  ] 


(19) 


and  the  phase  angle  by: 

0  *  arc  tan 


ct“ 


arc  tan 


ctco 

: — ? 

h°n 

'-O 


(20) 


in  which  the  fundamental  structural  resonant  frequency,  con,  is  given  by: 


Kt 

co  ■  2nf  -  (=-) 
n  n 


V2 


(21) 
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It  should  be  noted  that  this  frequency  cannot  be  determined  by  direct 
experimental  measurement  since  it  involves  the  Theodorsen  stiffness  and 
inertia  but  excludes  damping.  It  is  also  a  well  established  fact  that 
for  stationary  bodies  the  Strouhal  number 


V 


(22) 


is  a  unique  function  of  body  thickness  Reynolds  number,  Above  a 
moderate  value  of  *t(*  500  for  circular  cylinders)  the  Strouhal  number 
for  stationary  bodies  is  approximately  constant.  Thus,  with  the 
simplified  model  presented  above,  resonant  vibrations  are  expected  ft 
when  the  free  stream  velocity,  V,  approaches  the  value  given  by  V  *  -j— . 

t» 


It  has  been  discovered  from  previous  investigations  however  (1) 
that  the  Strouhal  number  falls  below  its  constant  value  in  the  vicinity 
of  structural  resonance  where  large  vibrational  amplitudes  prevail. 

This  fact  illustrates  the  influence  of  the  vibratory  motion  of  the  plate 
on  the  vortex  shedding  mechanism  and  thus  on  the  forcing  moment  coeffi¬ 
cient  nif0.  In  this  case  a  dependence  of  the  Strouhal  number  and  the 
moment  coefficient  on  the  motion  of  the  plate  can  lead  to  a  non-linear 
equation  of  motion.  The  linear  system  described  above  cannot  therefore 
be  expected  to  give  an  exact  representation  of  the  vibration  phenomenon 
but  it  can  be  used  as  a  first  approximation  for  the  evaluation  of  the 
forcing  moment  coefficient  from  experimental  data. 


V.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 
A.  System  Constants 

As  has  been  stated  in  Section  IV  the  coefficients  Is,  Ip,  cs, 

Cy,  and  Ks  are  independent  of  the  free  stream  velocity  and  time  and 
can  be  predetermined  and  are  summarized  in  Table  I. 

The  structural  moment  of  inertia,  Is,  for  any  spring-plate  system 
has  been  computed  from  the  geometry  of  the  system  (Table  I).  The  moment 
of  inertia  of  the  plate  itself  accounts  for  9/10  of  the  total  structural 
moment  of  inertia. 

In  order  to  determine  the  natural  frequencies  and  damping  coeffi¬ 
cients,  free  oscillation  tests  have  been  performed  in  air  (ct  ■  cs)  and 
in  still  water  (c3v  *  cs  +  Cy)  and  are  described  in  detail  in  Appendix  B. 
Such  tests  were  performed  for  each  plate-spring  system  at  the  beginning 
of  the  vibrational  tests  and  for  many  systems  at  the  end  of  the  vibrational 
tests  to  check  the  dependability  of  the  results.  The  natural  frequencies 
in  air  at  the  end  of  the  vibrational  tests  are  in  all  cases  slightly 
smaller  than  at  the  beginning  of  the  tests.  The  damping  in  air  at  the  end 
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of  the  vibrational  tests  i  -  considerably  higher  than  at  the  beginning. 

It  should  be  kept  in  mind  that  at  the  beginning  of  the  vibrational  tests 
the  whole  mounting  assembly  had  been  cleaned  and  the  bearings  lubri¬ 
cated.  At  the  end  of  the  vibrational  tests  the  lower  portion  of  the 
shaft-bearing  system  was  submerged  in  a  small  cavity  cf  water  and  some 
dirt  had  accumulated  over  uhe  bearings.  Thus,  apparent  moment  of 
inertia  effects  due  to  submergence  of  the  lower  shaft  can  account  for 
the  smaller  natural  frequencies  at  the  end  of  the  tests  and  the 
accumulated  dirt  around  the  bearings  for  the  increased  damping.  In 
all  cases  the  structural  damping  is  very  small;  about  l/o  of  the 
critical.  The  natural  frequency  in  still  water  remained  constant 
during  the  course  of  vibrational  tests.  A  small  difference  is  obtained 
only  in  the  plate  D  -  spring  II  system  possibly  due  to  a  loosening  of  c.sv 
the  spring-shaft  connection.  The  relative  total  damping  coefficient, 
in  still  water  is  small  b\it  the  absolute  viscous  damping  coefficient,  c 
Cy,  is  larger  than  the  corresponding  structural  coefficient,  cs. 


The  structural  spring  coefficient,  Ks,  was  obtained  from  in-place 
static  calibration  as  described  in  Appendix  B.  The  apparent  moment  of 
inertia.  Ip,  for  the  vibr sting  plate  can  be  predicted  theoretically  from 
potential  flow  theory  as  nas  been  seen  earlier  (Equation  8) .  This 
formula  gives  an  apparent  moment  of  inertia  for  the  tested  plates: 

Ip  -  36.5?  x  lO"1*  lb. in. sec2. 


Since  in  the  present  tests  other  structural  parts  of  the  plate-spring 
system  are  submerged,  the  true  apparent  moment  of  inertia  should  be 
larger  than  this  theoretical  value.  The  apparent  moment  of  inertia  of 
the  plate-spring  system  is  easily  computed  from  the  natural  frequencies, 
fn,  determined  experimentally  in  still  air  and  in  still  water.  Thus 
we  have: 


2nf 


n  air 


co 


n 


air 


s 


(23) 


% 


2nf 


n  still  water 


co 


n  still  water 


I  +1 
s  p 


(2h) 


and 


n  air 


n  still  water 


2 

■)  -  1) 


(23) 


The  values  of  apparent  moments  of  inertia  computed  in  this  fashion  are 
tabulated  for  each  plate-spring  system  in  Table  I.  They  agree  very 
closely  with  the  '-alue  predicted  theoretically,  being  in  general  slight¬ 
ly  larger. 
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B.  Exploratory  Tests 


1.  Experimental  Procedure 

Exploratory  tests  were  performed  with  plate  C  -  spring  I  to 
reveal  the  pertinent  characteristics  of  the  plate  motion  and  of  the  wake 
structure . 

The  ambient  free  stream  velocity  range  yielding  measurable 
vibrational  excitation  of  the  plate  was  first  defined.  The  following 
measurements  were  then  made  at  several  discrete  velocities  selected  so 
as  to  cover  this  range: 

a)  Vibrational  amplitudes  of  all  the  significant  harmonic 
components  of  the  plate  motion. 

b)  Mean  velocity  distribution  across  the  half  width  of  the 
wake  at  a  longitudinal  station  x  ■  l.$t  downstream  of  the  trailing  edge. 

c)  Fluctuating  total  head  spectral  distribution  at  station 
x  ■  1.5t  and  at  the  transverse  position  y/t  at  which  the  discrete  compo¬ 
nent  of  the  spectrum  was  a  maximum.  (See  Fig.  13  for  coordinate  reference). 
The  results  obtained  in  these  exploratory  tests  are  presented  and  discussed 
in  the  following  sections. 


Plate. 


2.  Vibrational  Characteristics 

The  velocity  range  of  significant  vibrational  excitation  is 
from  5  fps  to  9  fps  for  the  system  tested.  Throughout  the  whole  range 
of  excitation  the  motion  of  the  plate  observed  on  the  screen  of  the 
oscilloscope  had  the  appearance  of  a  sinusoid  with  continuously  varying 
amplitude.  For  five  velocities  within  the  range  of  excitation  the 
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frequencies,  fa,  of  vibration  and  their  corresponding  RMS  amplitudes,  a, 
were  picked  up  with  the  vibration  analyzer.  The  results  are  tabulated 
in  Table  II.  In  all  cases  two  frequencies  were  present,  the  fundamental, 
fao»  and  the  third  harmonic,  fa, .  The  RMS  amplitude  of  the  third  harmonic, 
03,  is  much  smaller  than  that  or  the  fundamental,  3q,  the  ratio  of  ampli¬ 
tudes,  03/00,  increasing  with  increasing  velocity  and  fundamental  frequency. 

3.  Wake  Characteristics 


Mean  velocity  wake  transverses  were  made  for  the  same  five 
ambient  velocities  with  the  plate  vibrating.  A  trial  was  made  to  dupli¬ 
cate  the  ambient  flow  conditions  (mean  velocity  and  mean  static  pressure) 
which  were  obtained  with  the  plate  clamped  stationary.  The  duplication 
of  flow  velocity  was  not  very  successful  due  to  operational  difficulties. 
Therefore  the  mean  velocity  traverses  for  the  stationary  plate  do  not 
correspond  to  exactly  the  same  mean  ambient  velocity  as  for  the  vibrational 
tests.  The  mean  velocity  distributions  obtained  are  presented  in  Fig.  Al. 
In  all  cases*  the  velocity  defect  near  the  centerline  of  the  wake  is  larger 
for  the  vibrating  than  for  the  stationary  plate.  Within  the  accuracy  of 
the  measurements,  the  width  of  the  wake  is  not  affected  by  the  motion  of 
the  plate,  being  in  general  the  same  in  both  cases. 

Spectral  analyses  of  the  total  head  fluctuations  were  made  for 
four  ambient  velocities  with  the  plate  vibrating.  These  spectra  were 
obtained  at  x  -  l.$t  and  at  the  y/t  at  which  the  fundamental  discrete 
component  had  maximum  energy  density.  The  normalized  spectral  distri¬ 
butions  are  presented  in  Figs.  A 2  through  A$.  These  spectra  illustrate 
the  turbulent  nature  of  the  early  wake  and  the  strong  discrete  periodici¬ 
ties  present.  The  frequencies  at  which  discrete  energy  is  concentrated 
in  the  spectra  are  tabulated  together  with  the  corresponding  vibrational 
frequencies  in  Table  II.  It  is  seen  that  in  all  spectra  a  fundamental 
discrete  frequency  is  present  along  with  its  second  and  third  harmonics. 

The  fundamental  frequency  is  exactly  the  same  as  the  fundamental  frequency 
of  plate  vibrations.  For  this  reason  we  will  henceforth  omit  the  sub¬ 
scripts  distinguishing  the  wake  and  vibrational  frequencies. 

An  estimate  of  the  discrete  energy  was  made  by  computing  the  area 
of  the  spectral  peak  at  the  fundamental  discrete  frequency  and  subtrac¬ 
ting  the  continuous  turbulent  background.  The  energy  in  the  fluctuating 
total  head  at  the  fundamental  discrete  frequency  can  be  seen  (Table  II) 
to  account  for  roughly  5Q^o  of  the  total  energy  of  the  spectrum  at  this 
location.  The  total  head  fluctuation  spectra  presented  here  for  the 
vibrating  plate  show  only  quantitative  differences  when  compared  with 
the  spectra  obtained  for  the  same  plate  clamped  stationary  in  the  flow 
(Ref.  1*) .  These  latter  spectra  are  also  shown  in  Fig.  A2  and  Fig.  Ai*  to 
facilitate  comparison. 
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TABLE  II  RESULTS  OF  EXPLORATORY  TESTS 
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U.  Conclusions 


The  conclusions  drawn  from  the  preliminary  tests  served  as  a 
guide  for  the  final  tests. 

a)  The  mounting  assembly  and  sensing  elements  for  recording 
the  motion  of  the  plate  proved  satisfactory  in  all  respects. 

b)  Since  the  amplitude  of  the  third  vibrational  harmonic  is 
substantially  smaller  than  the  amplitude  of  tne  fundamental,  the  motion 
of  the  plate  may  be  considered  simple  harmonic  for  the  purpose  of  further 
analysis . 


c)  The  RMS  of  the  vibrational  amplitude  represents  to  an 
extremely  good  degree  of  accuracy  the  RMS  of  the  amplitude  of  the  funda¬ 
mental  vibrational  frequency. 

d)  The  large  percentage  of  the  discrete  energy  in  the  total 
head  fluctuation  spectra  indicates  that  the  RMS  of  the  total  head  fluc¬ 
tuation  can  reflect  in  a  qualitative  manner  changes  in  the  discrete  energy. 

e)  The  similarity  of  the  spectral  measurements  for  the 
stationary  and  vibrating  cases  substantiate  the  concept  that  vibrations 
are  initiated  by  pressure  fluctuations  resulting  from  the  shedding  of 
vortices  in  the  wake. 

f )  The  fundamental  frequency  of  total  head  fluctuation  in 
the  wake  is  the  same  as  the  fundamental  frequency  of  plate  vibration. 

C.  Final  Tests  -  Vibrational  Measurements 


1.  Experimental  Procedure 

For  each  plate-spring  system  the  torsion  spring  was  cali¬ 
brated  in  place.  The  natural  frequency  as  well  as  the  damping  coeffi¬ 
cient,  cB,  of  the  system  in  still  air  (water  tunnel  empty)  and  in  still 
water,  csv  ■  cs  +  Cy,  were  determined  experimentally  with  free  oscilla¬ 
tion  tests.  The  ambient  free  stream  velocity  was  increased  by  small 
increments  AV  and  the  following  measurements  were  made. 

a)  The  mean  ambient  velocity  -using  the  calibrated  contraction 
meter  of  the  water  tunnel. 

b)  The  frequencies  of  the  vibrational  motion  of  the  plate 
using  the  vibration  analyzer, 

c)  The  RMS  of  the  vibratory  motion  of  the  plate  with  the 
AC  Balantine  voltmeter. 
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The  velocity  was  increased  up  to  the  point  at  which  measurable 
vibrations  ceased.  The  velocity  was  then  decreased  by  small  decrements 
AV  and  the  same  quantities  were  measured. 

The  above  measurements  were  performed  with  the  following 
plate-spring  systems. 

Plate  A  -  spring  I 

Plate  B  -  spring  I,  II 

Plate  C  -  spring  I,  II 

Plate  D  -  spring  I,  II 

For  plate  D  -  spring  I,  II  and  plate  B  -  spring  I,  the  natural  frequency 
and  damping  coefficient  in  air  and  in  still  water  were  redetermined  ex¬ 
perimentally  at  the  end  of  the  tests. 

In  all  cases  the  vibratory  motion  of  the  plate  was  observed  on 
the  screen  of  the  oscilloscope  and  the  qualitative  nature  of  the  motion 
was  recorded. 

2.  Qualitative  Vibrational  Characteristics 

The  fundamental  vibrational  frequency,  f0(cps),  as  given  by 
the  vibration  analyzer  and  the  RMS  of  the  vibrational  amplitude,  a^(rad), 
are  plotted  versus  the  ambient  velocity,  V(fps),  for  each  plate-spring 
system  in  Figs.  A6  through  A9. 

The  qualitative  character  of  the  vibrations  is  discussed  in  the 
present  section. 

Plate  B 


Measurable  vibrations  start  at  a  certain  velocity  with  a  fre¬ 
quency  smaller  than  the  fundamental  natural  frequency  of  the  plate¬ 
spring  system.  For  a  small  range  of  increasing  velocity  the  fundamen¬ 
tal  vibrational  frequency  varies  linearly  with  the  ambient  free-stream 
velocity  corresponding  to  a  constant  value  of  the  Strouhal  number, 

S-t  *  f0t/V  “0.258.  The  vibrational  amplitude  increases  with  increasing 
velocity  and  frequency  over  this  range.  In  this  range  the  only  vibra¬ 
tional  frequency  present  is  the  fundamental.  Only  in  the  upper  end  of 
this  range  does  an  extremely  small  (hardly  detectable)  third  harmonic 
make  its  appearance.  The  vibrational  frequency  in  this  range  is  constant 
(called  "stable"  hereafter)  as  far  as  could  be  determined  experimentally 
but  the  amplitude  varies  (unstable) .  The  variation  of  the  amplitude  was 
reflected  in  fluctuations  of  the  needle  of  the  voltmeter  and  was  observed 
also  on  the  screen  of  the  oscilloscope. 
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With  a  small  increase  of  velocity  above  this  range  the  amplitude 
of  vibrations  increases  sharply  and  the  fundamental  vibrational  fre¬ 
quency  falls  below  the  value  predicted  by  the  constant  Strouhal  number 
of  the  earlier  range.  A  considerable  range  of  velocities  exist  where 
large  amplitudes  are  obtained  while  the  Strouhal  number  decreases  con¬ 
tinuously  with  increasing  velocity.  The  third  harmonic  is  present 
throughout  this  range  with  increasing  amplitude  as  the  velocity  in¬ 
creases.  In  this  range  both  the  frequency  and  the  vibrational  amplitude 
are  stable  (see  Figs.  A10  and  All).  Audible  noise  was  present. 

Increasing  the  velocity  further,  a  velocity  is  reached  at  which 
the  vibrations  become  intermittent  with  amplitude  variable  as  much  as 
10  fold  (see  Figs.  A12  and  A13) .  The  amplitude  of  the  third  harmonic 
seems  to  attain  its  maximum  value  at  this  point. 

For  still  larger  velocities,  very  small  vibrations  are  obtained 
at  frequencies  indicating  an  increasing  Strouhal  number  with  increasing 
velocity.  Thus  three  velocity  ranges  can  be  distinguished,  as  was  re¬ 
ported  previously  by  Ippen  et  al.  (1) : 

(1)  Range  I  in  which  the  Strouhal  number  remains  constant 
with  sinusoidal  vibrations  of  varying  amplitude. 

(2)  Range  II  in  which  the  Strouhal  number  falls  below  its 
constant  value  with  large  stable  vibrational  amplitudes  and  increasing 
amplitude  of  the  third  harmonic  with  increasing  velocity. 

(3)  Range  III  in  which  small  varying  amplitude  vibrations  are 
obtained  and  increasing  Strouhal  numbers.  The  vibrational  characteris¬ 
tics  described  above  are  exhibited  by  plate  B  with  both  springs. 

Plate  D 


Plate  D  exhibits  vibrational  behavior  which  is  qualitatively 
similar  to  that  of  plate  B.  For  plate  D,  however,  there  is  a  more 
gradual  change  in  amplitude  as  the  velocity  is  increased  from  range  I 
to  range  II.  Vibrational  measurements  in  the  range  III  were  not  made 
with  plate  D.  The  behavior  of  plate  D  in  ranges  I  and  II  is  illustra¬ 
ted  in  Figs.  All)  through  A17. 

Plate  C 


Throughout  the  whole  range  of  excitation  the  vibrational  ampli¬ 
tude  was  unstable  (see  Figs .  A18  and  A19) .  The  third  harmonic  was 
present  but  hardly  detectable.  Due  to  the  small  vibrational  ampli¬ 
tudes  obtained  for  this  plate,  the  range  of  constant  Strouhal  number 
(Range  I)  could  not  be  measured. 
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Plate  A 


Extremely  small  vibrational  amplitudes  were  obtained.  Throughout 
the  whole  range  of  measurable  vibrational  excitation  the  amplitude  was 
unstable  and  the  Strouhal  number  decreased  slightly  with  increasing 
velocity. 

3.  The  Strouhal  Number. 

Each  plate  exhibits  its  own  variation  of  vibrational  fre¬ 
quency  with  ambient  free  stream  velocity,  indicating  the  importance  of 
trailing  edge  geometry  in  the  vortex  shedding  mechanism.  The  Strouhal 
number  was  computed  for  plates  B,  D,  tested  with  both  springs,  from  the 
data  taken  in  the  range  of  linear  dependence  of  the  fundamental  vibra¬ 
tional  frequency  on  the  velocity.  An  average  value  of  the  Strouhal 
number  was  defined  for  plate  A  from  the  data  taken  throughout  the  whole 
range  of  excitation.  The  Strouhal  number  was  not  computed  for  plate  C 
due  to  the  lack  of  experimental  data  in  the  velocity  range  I. 

The  Strouhal  numbers  for  the  three  plates  are  tabulated  in 
Table  III  together  with  the  values  obtained  for  the  same  plates  by 
other  investigators .  The  Strouhal  numbers  on  lines  (3) (U)  are  based 
on  the  fundamental  frequency  of  total  head  fluctuation  in  the  wake 
of  stationary  plates  at  a  velocity  of  about  10  ft/sec.  The  values  on 
line  (3)  agree  very  well  with  the  present  values,  the  difference  being 
within  the  experimental  error.  The  value  on  line  (1*)  differs  appreci¬ 
ably  from  that  obtained  for  the  present  tests.  This  difference  may 
result  from  the  rougher  plate  surface  and  lower  free-stream  turbulence 
level  in  the  present  tests.  The  values  in  line  (5)  based  on  funda¬ 
mental  vibrational  frequencies  are  not  in  good  agreement  with  the 
present  values  although  the  trend  is  the  same.  Differences  in  plate 
allignment  and  free-stream  turbulence  level  may  be  the  cause  of  the 
noticed  discrepancies. 

The  values  of  the  Strouhal  number  for  plates  A,  B,  and  D  differ 
substantially,  reinforcing  earlier  findings  (Ippen  et  al.  (1), 

Heskestad  (8))  that  the  trailing  edge  geometry  influences  the  vortex 
shedding  frequency  to  a  great  extent.  This  is  understandable  since  the 
time  required  for  the  vortices  to  mature  and  be  detached  as  discrete 
entities  from  the  separated  shearing  layers  depend  strongly  on  the 
early  wake  boundary  conditions.  The  work  of  Gongwer  (11)  and 
Eagleson  et  al.  (ii)  shows  the  importance  of  separation  point  boundary 
layer  thickness  on  the  Strouhal  number.  Thus,  until  the  detailed 
mechanism  of  vortex  formation  and  discharge  is  understood,  a  universal 
Strouhal  number  cannot  be  formulated. 
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TABLE  III 


Strouhal  Number 


Plate  A 

Plate  B 

Plate  C 

Plate  D 

1. 

Plate’  v 

E=2 

CD 

CD 

2. 

Present  tests 

0.258 

0.213 

h- 

Ref.  (3) 

0.237 

0.251 

0.223 

k_ 

Ref.  (2) 

0.255 

Ref.  (1) 

0-^9 _ 

0-27? _ 

0.21:6 

0.230 

It  has  been  stated  in  section  V-C-2  that  the  vibrational  Strouhal 
number  decreases  with  increasing  velocity  in  the  velocity  range  II  and 
then  increases  again  in  the  velocity  range  III.  The  maximum  and  minimum 
values  of  the  vibrational  Strouhal  number  obtained  for  the  different 
plates  are: 


Plate  A  S.  mnv  -  0.251 
t  max 

Plate  B  S.  „  -  0.258 
u  max 

Plate  C  S.  mjnr  -  0.257 
u  max 

Plate  D  S.  mav  -  0.213 
u  max 


St  min 
^t  min 
St  min 
^t  min 


0.231: 

0.185 

0.211: 

0.106 


The  larger  the  drop  in  Strouhal  number  from  range  I  to  range  II  the 
larger  the  range  of  velocities  in  which  large  vibrational  excitation 
occurs.  It  is  in  the  velocity  range  n  where  large  amplitudes  are  ob¬ 
tained.  It  is  thus  suggested  that  the  decrease  of  Strouhal  number  in 
the  range  II  is  the  result  of  dynamic  changes  imposed  by  the  motion  of 
the  plate.  It  is  the  plate-spring  system  rather  than  the  flow  velocity 
which  controls  the  frequency  in  this  range.  The  decrease  in  Strouhal 
number  seems  too  large  to  be  attributed  to  the  small  possible  trans¬ 
verse  shift  of  separation  points  caused  by  the  plate  motion. 


The  frequency  versus  velocity  relation  in  range  II  for  each 
plate  shows  a  similarity  for  the  two  tested  springs.  It  seems  to  be 
independent  of  the  absolute  magnitude  of  the  vibrational  amplitudes . 
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Of  Interest  is  that  for  plates  A,  B,  and  C  the  fundamental 
vibrational  frequency  obtained  at  a  certain  velocity  seems  to  be  inde¬ 
pendent  of  the  way  the  velocity  was  reached  (increasing  or  decreasing 
velocity) .  For  plate  D  an  hysteresis  effect  is  obvious  for  both  springs . 
With  decreasing  velocity  the  fundamental  vibrational  frequency  at  a 
certain  velocity  is  smaller  than  that  obtained  with  increasing  velocity. 
The  differences  are  small,  in  general  less  than  1  cps  and  are  more  pro¬ 
nounced  in  the  vicinity  of  the  lower  end*  of  range  II  while  diminishing 
on  the  upper  end.  It  should  be  noted  that  plates  A,  B,  and  C  have  more 
or  less  fixed  separation  points  while  the  trailing  edge  of  plate  D  gives 
the  possibility  of  forward  or  backward  shift  of  the  separation  points. 

It  is  thus  more  likely  that  the  frequency  hysteresis  effect  is  due  to 
possible  shift  of  the  separation  points. 

U.  Vibrational  Amplitude  versus  Ambient  Velocity 

In  Figs.  A6  through  A9  the  RMS  of  ohe  vibrational  amplitude , 
aJJ(rad),  is  plotted  versus  the  ambient  velocity,  V(fps).  From  these 
figures  the  effect  of  trailing  edge  geometry  on  the  vibrational  ampli¬ 
tudes  and  velocity  range  of  excitation  is  strikingly  clear.  Plates  B 
and  D  exhibit  maximum  amplitudes  of  the  same  order  of  magnitude  and 
much  larger  than  for  plates  C  and  A.  A  common  characteristic  of  plates 
B  and  D  is  that  the  trailing  edge  extends  into  the  early  wake  beyond 
the  separation  points. 

In  all  cases  the  maximum  amplitude  for  every  plate  is  smaller 
when  tested  with  spring  II  than  with  spring  I. 

The  amplitude  distributions  are  dissimilar  when  compared  from 
plate  to  plate  but  show  a  marked  similarity  for  the  same  plate  and  the 
two  springs  tested.  Thus  a  unique  presentation  of  the  results  in  a 
non-dimensional  form  is  possible.  Figs.  A20  through  A22  present  a 
single  curve  for  each  plate  of  the  dimensionless  amplitude,  7T0/ao  mp„, 
versus  the  dimensionless  velocity,  V/Vr,  Vr  being  the  velocity  at  wMch 
o0  max  occurs  for  each  plate-spring  system.  This  presentation  suggests 
the  unique  importance  of  the  velocity  and  trailing  edge  geometry  on  the 
vibrational  behavior  of  each  plate. 

A  careful  examination  of  the  results  obtained  indicates  that  the 
vibrational  amplitude  of  each  plate-spring  system  at  a  certain  velocity 
is  practically  independent  of  the  way  the  velocity  was  reached.  Some 
differences  may  exist  at  the  transition  from  range  I  to  range  II. 

<y.  Reproducibility  of  Results 

Confidence  in  the  quantitative  value  of  experimental  measure¬ 
ments  is  lost  if  the  measurements  are  not  repeatable.  Thus  upon  com¬ 
pletion  of  the  final  experiments  the  plate  B  -  spring  II  system  was  re¬ 
tested.  Natural  frequencies  and  damping  coefficients  for  the  two  tests 
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with  this  system  have  already  been  presented  in  Table  I  with  a  difference 
less  than  the  estimated  experimental  error.  The  vibrational  test  results 
are  presented  in  Fig.  A23  together  with  the  results  obtained  previously. 

No  comment  is  necessary. 

6.  The  Response  Curves 

In  Figs.  A2li  through  A27  the  vibrational  results  are  presented 
in  the  more  conventional  form  of  amplitude  response  curves.  The  RMS  of 
the  vibrational  amplitudes,  a0(rad),  are  plotted  versus  the  dimensionless 
frequency,  fQ/f_.  The  maximum  amplitude  for  all  plate-spring  systems 
occurs  in  the  vicinity  of  the  corresponding  natural  frequency,  in  general 
at  a  frequency  slightly  larger  than  the  natural  one.  The  figures  show 
a  rapid  rise  of  amplitude  left  of  resonance  for  plates  B  and  C  which  is 
especially  noticeable  with  spring  I.  No  such  behavior  is  shown  for  plate 
D  for  which  the  increase  of  amplitude  is  more  gradual.  For  all  plates 
the  amplitudes  obtained  in  the  vicinity  of  and  to  the  "right"  (higher 
frequencies)  of  resonance  are  independent  of  the  way  the  particular  velo¬ 
city  was  reached.  It  should  be  noted  that  stable  vibrations  were  ob¬ 
tained  in  this  region  for  plates  B  and  D.  At  the  "left"  (lower  fre¬ 
quencies)  of  resonance  larger  amplitudes  are  obtained  in  general  with 
decreasing  velocity  than  with  increasing  velocity.  This  is  especially 
noticeable  with  plate  D  for  which  the  frequency  also  exhibited  an 
hysteresis  effect.  This  region  corresponds  to  velocity  range  I  and  the 
transition  from  range  I  to  II  where  varying  amplitude  vibrations  were 
obtained. 


In  order  to  eliminate  the  explicit  role  of  velocity  in  the 
forcing  moment  expressed  as  M0  -  mf0  1/2  pH(2b)2V2,  the  experimental 
data  are  presented  in  Figs.  A28  through  A31  in  a  reduced  dimensionless 


form  of  C5o/72)/(30/72)maX  vs.  fD/fn.  T*1®  reduced  response  curves  in¬ 
dicated  that  for  all  plates  the  reduced  maximum  amplitude  occurs  at  the 
natural  frequency  of  the  system,  small  deviations  being  within  the  ex¬ 
perimental  error.  The  general  characteristics  of  the  reduced  response 
curves  are  otherwise  the  same  as  those  presented  in  Figs.  A2k  through 
A27.  The  hysteresis  effect  is  now  more  clearly  shown. 


7.  The  Forcing  Moment  Coefficient 


Equation  (19)  served  as  the  basis  for  the  computation  of 
the  forcing  moment  coefficient  mf0  from  the  obtained  experimental  re¬ 
sults. 


mfQ  i  pH(2b)2V2 


*tV 


2s2.  2  2-,-te 

■«o  )  +ct  coo  ] 


(19) 
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2  *8  1  2  2 

in  which  ■  (-| — E) ,  K  *  ^npb  HV  ,  and  to0  is  the  fundamental 

vibrational  frequency.  An  evaluation  of  the  term  Kp  for  the  range  of 
velocities  producing  measurable  vibrational  excitation  in  the  present 
tests  indicates  that  Kp«Xs.  Thus,  Kp  can  be  neglected  and  O' K_/I^, 
being  the  undamped  natural  frequency  of  the  plate-spring  system  in 
still  water.  The  error  introduced  in  Un  by  neglecting  the  Kp  term  is 
of  the  same  order  as  the  experimental  error  in  determining  tne  natural 
frequency  in  still  water  (i.e.  ±  1  cps) .  The  total  damping  coefficient 
is  given  by: 

ct  "  csv  +  I  nPb3HV •  (18b) 

csv  is  the  structural  plus  viscous  damping  coefficient  of  the  system  and 
was  assumed  independent  of  the  velocity  V  and  the  motion  of  the  plate. 

It  is  assumed  equal  to  the  damping  coefficient  obtained  by  the  free 
oscillation  tests  in  still  water  (see  Table  I).  The  term (9/i^npb3HV  re¬ 
presents  the  damping  coefficient  introduced  in  the  system  by  the 
Theodorsen  solution  of  the  flat  plate  oscillating  in  a  uniform  stream 
of  incompressible  fluid.  Evaluation  of  this  term  indicates  that  it  is 
of  great  importance  and  is  much  larger  than  the  term  csv. 

Letting  Cp  -  (9/hVipb^HV  ■  CV, 
we  can  write:  o  o  o  2  2  o 

ao[It  (“n  -“o  }  +  (csv+CV)  ®o  1 

mf°  "  — 

i  pH (2b) V 

All  the  quantities  in  the  right  member  of  the  above  expression  are  known 
and  mf0  can  be  easily  computed^  For  each  plate-spring  system  a  mean 
smooth  curve  was  drawn  in  the  a0  versus  fD  and  f0  versus  V  graphs  and 
from  these  curves  values  of  and  V  were  obtained  for  use  in  the  compu¬ 
tation  of  mf0  at  selected  frequencies.  (Use  of  a  smoothed  curve  rather 
than  the  actual  data  points  is  a  convenience  which  will  be  justified 
later).  Since  the  RMS  of  a0  (i.e.  a0  *  a0/\/2~ )  was  readily  available, 
the  RMS  of  mf0  (i.e.  mfp  *  mf0/  \J2~  )  was  computed.  Calculations  of 
If0  were  performed  at  equal  intervals  of  A(<on/co0)  at  the  left  of  resonance 
and  of  A  (coo/con)  at  the  right  of  resonance. 

Computations  of  nu> 0  were  made  for  plates  B,  C,  D,  and  both 
springs.  The  experimental  results  for  plate  A  were  not  considered 
accurate  enough  for  quantitative  correlations  since  in  this  case  the 
electrical  signal  representing  a0  was  of  the  same  order  of  magnitude  as 
the  background  noise  of  the  electronic  instruments  used  for  vibrational 
measurements.  For  every  frequency,  co0,  at  which  mf0  was  computed,  the 
corresponding  quantity  2bco0S0/V  was  also  computed.  This  dimensionless 
quantity  represents  the  ratio  of  the  RMS  transverse  velocity  of  the  plate 


(26) 

V2 

—  (27) 
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At  the  trailing  edge  to  the  ambient  free  stream  velocity.  For  all  plate- 
spring  systems  the  computed  m^0  is  plotted  versus  the  corresponding 
V*  *  2ba0a0/V  quantity  and  is  shown  in  figures  A32,  A33. 


The  correlation  of  mf0  with  2bea0a0/T  indicates  characteristics 
of  great  importance  for  the  understanding  of  the  vibrational  problem. 

A  large  increase  of  mf0  occurs  with  increasing  values  of  2hco030/7  indi¬ 
cating  a  strong  interaction  between  forcing  moment  and  motion  of  the 
plate.  For  small  values  of  2bco030/7>  at  left  of  resonance  and  in  the 
velocity  range  I  where  the  Strouhal  number  remains  constant  for  plates 
B,  D,  mf0  varies  linearly  with  2bco0a0/V .  Thus  mf0  in  this  range  can 
be  written: 


(■ 


2bco  a 
o  On 

V 


(28) 


The  values  of  and  are  essentially  the  same  for  both  plates  B,  D, 
with  either  spring.  The  value  of  is  of  great  physical  importance  in¬ 
dicating  that  with  no  vibrations  small  forcing  moments  are  present  which 
initiate  the  motion  of  the  plate.  In  the  vicinity  of  resonance  and  to  the 
right  of  resonance,  in  the  velocity  range  II  where  the  Strouhal  number 
falls  below  its  initial  constant  value,  the  rate  of  Increase  of  mf0  with 
increasing  2bayi0/V  is  decreased.  Here  too,  the  mfD  show  a  linear  depen¬ 
dence  on  2bca0ao/7  but  with  a  smaller  slope  than  in  range  I.  In  this  range 
the  nif0  can  be  written: 


"II 


2bco  a 

+  C?  (——> 

2n  7 


(29) 


For  both  plates  B,  D  and  for  both  springs,  the  C-j^  and  C2jj  values  are 

essentially  the  same.  The  experimental  results  indicate  a  larger  scatter 
for  plate  B  than  for  plate  D.  This  is  mainly  due  to  the  sharper  response 
curve  of  plate  B  and  to  larger  error  introduced  in  the  computation  of 
mfo  from  the  experimental  inaccuracy  of  cojj.  For  plate  C  no  experimental 
results  were  available  in  range  I.  In  range  II  the  slope,  is  exact¬ 

ly  the  same  as  that  for  plates  B  and  D  while  the  value  C]_TT  indicates  a 
difference  between  the  two  springs. 


Disregarding  the  minor  differences  in  the  behavior  of  the  various 
plates,  it  can  be  said  that  two  ranges  of  mj£  can  be  distinguished 
corresponding  to  the  two  velocity  ranges,  I,  II,  of  vibrational  excitation: 


2bco  a 
,  o  o 


) 


Range  I: 


V 


(28) 


f 


2boo i  o 

Range  II:  i—f*)  (29) 

with  C^j,  C2l>  and  C2  the  same  for  all  plates  and  springs;  and  C^j 
almost  the  same.  These  values  are  tabulated  in  Table  IV. 


we  have: 


Introducing  the  above  expressions  in  the  amplitude  relation  (19) 


Range  I:  a 


(C1  +C2  (~~-2))  f  «(2b)2  V2 
,  I  I  v  _ 

lIt2(“n2-“o2)^<cs,*CT)2“o2^ 


(30) 


Range  II:  a  * 


2bco  a 

(C  +C  (_^_£))  £  H(2b)^  V 
AII  II  V 


2  it  2 


2/..  2 


VV-V  +Kv  +  cv)  «0  3 


(31) 


TABLE  IV 

NOMINAL  VALUES  OF  MOMENT  COEFFICIENT  CONSTANTS 


Plate  Spring 

% 

% 

C1 

XII 

C2 

*11 

CI 

x  2 
lbs-sec 

CII 

lbs-sec 

B 

I 

0.0025 

2.425 

0.0090 

1.765 

O.OO783 

0.00569 

B 

II 

0.0025 

2.425 

0.0090 

1.765 

O.OO783 

0.00569 

C 

I 

0.0078 

1.765 

O.OO783 

0.00569 

C 

II 

— 

0.0042 

1.765 

O.OO783 

0.00569 

D 

I 

0.0025 

2.425 

0.0090 

1.765 

0.00783 

0.00569 

D 

II 

0.0025 

2.425 

0.0090 

1.765 

0.00783 

0.00569 

Theodorsen  Damping  Coefficient,  C  -  £  n  p  b%  ■  0.00569  lb. sec2  (Eq.  26) 


Solving  explicitly  for  aQ  we  have: 


2  „2 


Range  I: 


0.  §  H(2b)^  V 

XI  * 


(32) 
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(33) 


2  „2 


Range  II:  a  * 


Let: 


C2  |  H(2b)2  V' 


[It2(»n2-«o2)‘  *  (o„*  CV)'«02]'72-  C2n  |  H(2b)^Vo»0 


2.. 


and 


Then: 


£II 


H(2b)3  -  Cj 

(3U) 

H(2b)3  -  Cn 

(35) 

H(2b)2  V2 

Range  I:  a  «  -■ 

0  2/2  2s2 


[V(«bC—o‘>'*  -  °I  T“o 


(36) 


C2  |  H(2b)2  V2 

Range  II:  7T - 2 - - - tt -  (37) 

[It2(“n2-“o2>  ♦«.„*">  O*  ‘  CII  T“o 

From  these  two  expressions  it  is  obvious  that  the  coefficients,  C,  Cj, 
and  Cjj  have  the  same  dimensions.  A  quantitative  evaluation  of  C,  Cj, 
and  Cxi  yields  the  important  result  that  Cj  >  C  and  Cji  s  C.  (See 
Table  IV.)  The  two  last  expressions  for  a£  (equations  36  and  37)  will 
constitute  the  basis  for  further  discussion  of  the  vibrational  behavior. 

In  order  to  check  the  hysteresis  effect  on  the  relationship  of 
moment  coefficient  versus  2boi050/V,  the  following  computations  were 
made  for  plate  D-spring  II  for  which  the  hysteresis  is  more  pronounced. 

A  smooth  curve  was  drawn  in  the  RJ  versus  f0  graph  through  points  ob¬ 
tained  with  increasing  velocity,  and  a  second  one  through  points  ob¬ 
tained  with  decreasing  velocity.  The  mf0  and  the  corresponding  2bco0a0/V 
were  then  computed  using  points  on  the  two  smoothed  curves  and  they  are 
presented  in  Fig.  A3l*  together  with  the  results  obtained  from  the  mean 
curve.  No  substantial  difference  is  indicated  and  the  approach  followed 
in  the  mf£  computation  is  thus  justified. 

8.  Range  I  -  The  Frequency  Range  of  Increasing  Vibrational 
Amplitudes . 

It  has  been  derived  in  the  previous  section  that  in  range  I 
the  amplitude  is  given  by  the  expression: 


-36- 


a  * 

o 


(36) 


|  H(2b)2  V2 

/  .2 

+  (csv  +  CV)  to/]2 


CI  V“o 


In  this  range  the  Strouhal  number  remains  constant  for  each  plate.  Thus, 
the  velocity  is  given  by  the  expression: 


tot 

2^  * 


(38) 


Substituting  for  V  in  (36)  we  get: 


T 

o 


\  2  H(2b)2  (2rS~) 


[lt2(“n2  "“o2)2+  (csv  +  C(2^T))iifflo2]^  - 


(39) 


Let: 


and 


tot2  \ 

[It2  (“n2-“o2)2+(csv  +  C(2^))  O  "FlK} 


C_  (r4r-)  to  2  *  (to  ) 
I  v2nS/  0  2  K  o' 


0*0) 

0*1) 


Thus,  the  denominator  of  the  equation  (39)  is  (co0)  -  F2  (w0) .  Physical 
reality  makes  it  necessary  that  the  denominator,  F]_  (coG)-  F2  (too),  be 
larger  than  zero  so  that  finite  amplitudes  are  obtained.  Since  Cj  >  C, 
in  the  vicinity  of  resonance  F]_  (»0)  can  be  smaller  than  F2  (too)  depen¬ 
ding  on  the  value  of  csv.  In  this  case  the  relation  (36)  no  longer 
holds  and  some  dynamic  adjustment  must  occur.  Experiment  has  shown 
this  adjustment  to  involve  a  change  from  Cj  >  C  to  Cu  ■  C,  while  the 
Strouhal  number  falls  below  its  constant  value.  Thus,  we  have  a  classi¬ 
cal  Eigenvalue  problem  defined  by  the  relation: 


F1  (toQ)  -  ?2  («o)  -  0  (42) 

A  qualitative  representation  of  this  equation  is  shown  in  Figure  ll*» 

Depending  on  the  actual  value  of  csv  the  curve  Fi  (too)  may  inter¬ 
sect  the  curve  F2  (too)  at  two  points,  ©l  and  ©u  (Eigenvalues),  at  one 
point,  or  may  not  intersect  F2(©q)  at  any  point  for  large  values  of  csv. 
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i 


The  equation  F^  (aQ)  -  F^  (aQ)  -  0  is  written! 

[I.2(«  2 -a  2)2+(c  +C(^-))2»  2]^2  -  C,(*4-)«  2 

1  t  '  n  o  '  '  sv  ^2nS^"  o  1  Iv2nS^/  o 


Fig.  lii  Graphical  Representation  of  Equation,  F^(mq)  -F^i® Qh» 


0*3) 
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In  the  present  tests,  the  value  of  csv  is  much  smaller  than  the 
Theodorsen  term,  CV,  in  range  I  and  as  a  first  approximation  can  be 
neglected  for  the  general  discussion  of  equation  (1*3).  Thus, 
equation  (U3)  can  be  written: 


,  2/2  2 ^  ,  „2,  t  ^  1*  a»r  2/  U 

h  (“n  "®o  }  +C  (55Sr}  “o  CI  “o 


or,  solving  for  the  ratio 


co  ' 
(— ) 


CO. 


m 


(1*5) 


Since  Cj  >  C  two  real  solutions  for  coQ  are  obtained.  These  two 
Eigenvalues  are  called  the  lower  and  upper  critical  frequencies ,  co^  and 
coy,  and  hence  define  the  frequency  limits  of  range  II.  Using  the  above 
expression  and  the  actual  values  of  St  and  It  for  plates  B  and  D,  the 
col  a^  o>u  were  computed  for  both  springs .  The  values  of  coq  thus  obtained 
agree  very  well  (see  Table  V)  with  the  experimentally  determined  fre¬ 
quencies  which  bound  range  II.  These  experimental  critical  frequencies 
were  determined  from  Figs.  A7  and  A9.  The  lower  critical  frequency  and 
velocity  were  defined  by  the  point  at  which  the  Strouhal  number  first 
leaves  its  constant  value.  The  upper  critical  point  was  selected  by 
the  velocity  at  which  the  vibrational  amplitude  abruptly  decreased. 

While  the  relationship  7  *  co0t/2nSt  does  not  hold  for  co^,  the  value  ob¬ 
tained  from  the  above  equation  agrees  extremely  well  with  the  experi¬ 
mental  value.  Eq.  (U5)  thus  predicts  accurately  the  range  of  frequency 
over  which  large  amplitudes  are  obtained. 


It  is  of  great  interest  to  notice  here  that  the  importance  of 
the  trailing  edge  configuration  in  the  vibrational  behavior  enters  into 
equation  (1*5)  through  the  "steady-state"  Strouhal  number  St  and  Cj.  The 
larger  the  Strouhal  number  the  smaller  the  range  of  co„-col  where  large 
vibrations  occur.  The  nominal  value  of  Cj  is  independent  of  the  trailing 
edge,  however  wake  studies  (It)  have  implied  a  variability  in  the  steady 
state  moment  coefficients. 


An  exact  solution  of  equation  (1*3)  by  trial  incorporating  the 
actual  values  of  csv  for  each  plate-spring  system,  gives  values  of  col 
and  co^  only  slightly  different  from  the  ones  obtained  by  the  simplified 
equation  (1*5)  •  See  Table  V. 


Experimental  results  obtained  with  the  same  plates  reported  in 
Ref.  (1)  were  available  covering  a  much  larger  range  of  co^,  Using 
equation  (1*5)  and  the  values  of  «o_  and  from  Ref.  (1)  for  plates  B 

and  D  and  the  value  of  Cj  obtained  from  the  present  analysis,  values 
of  col  and  ®u  w#re  computed  and  compared  with  the  experimental  values  of 
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TABLE  V  CRITICAL  FREQUENCIES 


Equation  (1*5) _ Equation  (1*3) _ Experimental 


Plate 

Spring 

ESS 

fL  cps 

fu  cps 

B 

I 

66.9 

78.5 

67.6 

77.1 

69 

76.5 

B 

II 

69.3 

81.1 

70.3 

79.8 

7h 

80.8 

D 

I 

65 

79 

66 

77.2 

65.6 

79.2 

D 

II 

69 

83.9 

70.1* 

81.8 

69.1 

83 

frequencies  at  which  the  Strouhal  number  falls  below  its  constant  value 
left  of  resonance  and  the  point  where  vibrations  become  intermittent  at 
the  right  of  resonance.  All  the  results  are  summarized  in  Fig.  A35.  The 
predicted  and  experimental  values  are  in  extremely  good  agreement.  This 
agreement  suggests  that  the  value  of  Cj  obtained  from  the  present  tests 
performed  over  a  quite  narrow  range  of  natural  frequency  (f_  from  70  cps 
to  75  cps)  is  applicable  for  a  much  wider  range  of  natural  frequencies 
and  thus  of  velocities,  at  least  for  the  same  test  plates. 

9.  The  Critical  Value  of  csv 

It  has  been  discussed  in  the  previous  section  that  for  large 
values  of  csv,  the  curve  F2(co0)  may  not  intersect  the  curve  F^(co0) . 

Since  the  range  (II)  of  large  amplitude  vibrations  would  then  be  absent, 
it  is  of  practical  importance  to  determine  the  limiting  value  of  cBV  for 
which  the  equation  Fi(co0)  -F2(ea0)  ■0  has  only  one  real,  positive  root. 
Thi3  is  done  most  easily  by  returning  to  the  equation  of  motion,  Eq.  (15) 
and  substituting  the  assumed  solution,  Eq.  (17) .  Again  neglecting  0  (I) 
as  small  we  obtain: 


-Wv  ict"oV  KtV  t°li*c2i  -T-2'  SpbW.W 

Taking  the  imaginary  part  of  Eq.  (1*6)  and  solving  for  : 

2CiTpb2HV2  sin  0 

a - l -  (1*7) 

co  [c+  -CTVsin0] 

O  v  X 

The  denominator  of  Eq.  (1*7)  will  have  a  single  root  when  sin0-  1 

From  Eq.  (20)  this  can  be  seen  to  occur  when  co  ■  a  .  Thus,  the  minimum 

on 


I 


£ 

I 

I 


I 


!' 


value  of  CgV  in  order  to  avoid  large  amplitude  vibrations  is  obtained 
from: 

ct  -  CjV  (1*8) 


or 


csvc  "  ~  C) V 


(1*9) 


and  finally: 

C»VC  •  ^  <CI  -  c>  (5°) 

Referring  to  the  critical  damping  of  the  system,  cc  *21,  o^,  we 

have 


csvc  “  C) 

~  =  ^n  It  St 


(5D 


The  values  of  cSVc/cc  for  plates  B  and  D  have  been  computed 
using  Eq.  (51)  and  the  experimental  values  of  1^,  S-t  and  Ct.  They  are 
tabulated  in  Table  VI  in  comparison  with  the  relative  damping  actually 
experienced  in  the  tests. 


TABLE  VI  CRITICAL  VALUES  OF  csv 


csv/ cc 

CSVC/C| 

Plate  B  -  Spring  I 

O.82?o 

3.07?o 

Plate  B  -  Spring  II 

0.99?o 

3.00?o 

Plate  D  -  Spring  I 

1.50?o 

3.80?o 

Plate  D  -  Spring  II 

1.87?o 

3.90?o 

For  values  of  csv  >  cSVc  the  high  rate  of  increase  of  moment  co¬ 
efficient  with  the  self -excitation  parameter,  applicable  in  Range  I,  should 
persist  over  the  entire  range  of  vibrational  excitation.  The  Strouhal 
number  should  thus  retain  its  "steady-state"  value  and  due  to  the  higher 
damping,  much  smaller  vibrational  amplitudes  should  be  obtained. 

10.  Range  II  -  The  Velocity  Range  of  Large  Vibrational  Amplitudes 

The  development  of  section  V-C-8  has  made  possible  the  pre¬ 
diction  of  the  range  of  frequency  ©u_aL,  around  the  natural  frequency 
of  the  plate-spring  system,  for  which  large  vibrations  occur,  i.e. 
range  II.  To  complete  the  analysis,  the  functional  relationship  of  »0 
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and  V  within  this  range  of  frequency  must  be  established.  We  have  al¬ 
ready  seen  that  at  some  critical  value  of  the  free  stream  velocity,  the 
vibrational  amplitude  tends  to  infinity.  This  can  lead  to  one  of  two 
results : 


a)  structural  failure  resulting  from  unlimited  growth  of 
vibrational  amplitude. 

b)  large  but  finite  vibrational  amplitudes  resulting  from 
the  coming  into  play  of  some  amplitude  limiting  pheno¬ 
menon. 


The  present  experimental  results  indicate  clearly  that  the  ampli¬ 
tudes  in  Range  II  are  self -limiting,  which  implies  the  existence  of  some 
non-linear  damping  mechanisms  not  accounted  for  in  the  analytical  model 
in  use  here.  Nevertheless,  certain  aspects  of  the  Range  II  behavior 
which  have  considerable  engineering  significance  may  be  determined  with¬ 
out  a  detailed  understanding  of  the  amplitude-limiting  mechanism. 


From  vortex  street  considerations  Roshko  (hi*)  writes  the  mean 
body  circulation,  T^,  as 

*  *  y  ^ 


in  which  Vs  is  the  velocity  along  the  separation  streamline.  Assuming 
application  of  the  lifting  force  at  the  quarter-chord  point,  the  vortex- 
induced  moment  resulting  from  this  circulation  may  be  written 

V  2 

Hv  -  p  V  Tb  H  b/2  -  p  H  (g)  (^)  «  ifo  2pb2HV2  (53) 

or  - 

V  2 

("yO  “  St  ™fO 

When  the  plate  vibrates  in  response  to  this  vortex-induced 
forcing  moment,  the  velocity,  Vs,  along  the  separation  stream  line  will 
be  augmented  due  to  the  transverse  velocity  of  the  trailing  edge,  V<pg. 
For  bodies  of  constant  thickness,  Vs  should  be  of  the  general  form 

V!'’v2*vte 2  <5S> 

For  vibration  about  the  leading  edge  (as  is  the  case  being  considered 
here)  Vjg  is  given  by 

VTE  ■  (56) 
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Under  these  conditions 


V  2  2bco  a 

(-*)<-!+  — TT”  +  V* 


(57) 


When  transverse  vibration  does  not  exist  (V*  ■  0)  we  may  reasonably  ex¬ 
pect  (7S/V)Z  to  be  constant  and  Eq,  (5ii)  becomes,  in  Range  Is 


(f )  -  s  (t)  cii 

V  ax-0  1  t  1 


(58) 


Extrapolating  the  Range  II  behavior  to  V  ■  0  we  obtain  in  a  similar 
fashion: 

(f>  ■  st' (f)  Cln  m 

aII  U 


in  which 


St'  "  2n  TV  -  Vl7 

and  is  given  by  the  constant  Strouhal  number  of  Range  I  as 

V  ."fcL 

L  2nS, 


(60) 


(61) 


Under  the  assumption  that  the  proportionality  of  Eq,  (57)  is  constant 
at  least  in  the  absence  of  vibration: 

2 


V 

(-) 

V 


Oj  ■  0 


7 

(-) 

v 


(62) 


aH -0 


and  Eqs.  (58),  (59)  give: 


^lj 

st'  "  st  c r 


(63) 


II 


The  range  of  applicability 

“o  “  “L 
and  V  -  VL 

and  at  the  other  end  by 

co  ■  CO 

o  u 
and  V  ■  Vu 


of  Eq.  (63) 
(Predicted 
(Predicted 

(Predicted 


is  bounded  at  one  end  by: 
by  Eq.  U3) 
by  Eq.  6l) 

by  Eq.  h3) 


-1*3- 


Putting  eo0  ■  «#u  in  St',  Eq.  (63)  gives  a  predictor  for  Vu: 


V  ■  •• 

*  .1 


ft,*  r  -1  .  ClII 


u 


2^ 


[1  + 


CO 

’ 1)] 


m 


The  values  of  Vu  given  by  Eq.  (61;)  are  compared  below  with  the 
values  indicated  experimentally  for  plate3  5  and  0.  Vu  was  defined  ex¬ 
perimentally  as  the  velocity  at  which  the  vibrational  amplitude  abruptly 
decreases  (see  Figs.  A7  and  A9).  Note  that  when  the  nominal  values  of 
and  C]_jj  are  used  in  Eq.  (61*)  the  prediction  for  plate  D  is  in  con¬ 
siderable  error.  Examination  of  Fig.  A33  shows  that  the  actual  value  of 
Cij  for  plate  D  may  be  only  1/2  this  nominal  value  or  0.00125.  Predictions 
using  this  latter  value  are  also  ^iven  below. 


UPPER  CRITICAL  VELOCITY,  Vu 


Plate 

Spring 

Experimental  Vu 
.  fPs*  _ 

Theoretical  Vu 
(Eq.  61*)  fps. 

^1 

°in 

1. 

B 

I 

8.6 

8.2 

0.0025 

0.0090 

2. 

B 

II 

9.1 

8.5 

0.0025 

0.0090 

3. 

D 

I 

15.8 

10.3 

0.0025 

0.0090 

i*. 

D 

I 

15.8 

11*. 3 

0.00125 

0.0090 

5. 

D 

n 

16.6 

11.3 

0.0025 

0.0090 

6. 

D 

II 

16.6 

15.6 

0.00125 

0.0090 

The  predicted  frequency  -  velocity  relation  for  Range  II  using 
Eqs.  1*3,  6l,  63  and  61;  and  the  C's  of  lines  1,  2,  1;  and  6  above  are 
plotted  as  solid  lines  in  Figs.  A7  and  A9. 

Due  to  the  extreme  sensitivity  of  the  Range  II  behavior  to  the 
values  of  these  "steady-state"  moment  coefficients  and  to  avoid  the 
necessity  for  determining  them^an  alternative,  empirical  representation 
is  offered: 

It  was  observed  that  for  each  plate-spring  system,  the  value  of 
a0/V2  is  approximately  the  same  at  the  experimental  values  of  and  c^. 
This  is  shown  in  Table  VII. 
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TABLE  VII  CRITICAL  VELOCITY  CRITERION 


aQ/V2  at 


Plate  Spring  10^  rad  (sec2 /ft2) 


a0/V2  at  «u 
10^  rad  ( sec2/ft 2) 


B 

I 

6.8 

8.  it 

B 

II 

8  .it 

7.6 

D 

I 

2.75 

2.75 

D 

II 

2.0 

2.0 

Thus: 

=o  at  °L 

1  »<»>>  V 

[It2(“n2  -  “L2^  +  (csv  +  CVL)2  “L2^  * 

C  e  H(2b)2V  2 

■LII  £  u 

cc  at  o) 

0  u 

(I.  2(eo  2 -«  2)2  +  (c  +  CV  )2  to  2]^2  -  CV  co 

1  t  v  n  u  ’  v  sv  u'  u  J  u  u 

In  order 

that 

a„  at  co.  cE  at  co 

(-V)  ■  (■  -1--) 

\  Vu 

we  should  have: 

[I 

t2<“n2-“i,2)2*<0sv*CTL>2  “l2]V2-ctA 

* 

l  VV  -  “u2)2  *  (c»v  *  CT/  “»2)^  -  W»“U 

(66) 


(67) 


(68) 


With  predetermined  values  of  cor ,  co^,  and  V.,  the  above  equation  can 
be  solved  for  V„  since  the  left  hand  member  is  a  known  and  easily 
computed  quantity. 


Let 


tIt2<“n2-“L2)**<°»v*CV2  '  A 


(69) 


then 


*  ■  Ut2(»„2  -  %2)~  *  (0ST  ♦  CTU)-  <]  -  -  CVA 


2  2.% 


(70) 


or 


VV-02  * 

A  «  (c  +  CV  )<»  [1  +  -2 - - - S - ]  -  CV  to 

8V  UU  (c  +  CV)2£o2  UU 

'  sv  u  u 


(71) 


xtV-0 


For  values  of  —  ■  1  "  a  ^  «  1,  as  is  the  case  in  these  experiments, 

i  -  ^  mi  \  c 


(c  +CV  )  ® 

'  sv  u  u 


ItW>  ^  x  !t2(“n2  _“u2) 

[l+_t_n - u  3  »  i  +1  t  n  u 


(c  +  CV  )2<u  2 
'  sv  u'  u 


2  (c  +CV  )2£D  2 

'  ST  u'  U 


(72) 


or 


1 

A  "  W<Vu  +  I  -  C\% 

(c  +  CV  )  £0 
sv  u  u 


(73) 


and  finally 


po  9  c 

„  elrl  h  (“n  -“u  }  .  , 

V  -  ^  [  - - - - c ■  ) 

(A-  C  £0  )  £0 
sv  u  u 


(7U) 


The  above  expression  indicates  clearly  that  Vu  strongly  depends  on  csv. 
Thus  inaccuracies  in  determining  the  CsV  value  can  greatly  influence  the 
computed  value  of  Vy.  Using  equation  (7U),  the  values  of  0^  and  £0u  pre¬ 
dicted  by  equation  (L3)  and  the  experimental  values  of  It,  »n,  and  csv, 
the  Vu  vas  computed  for  plates  B  and  D  tested  with  both  springs.  Neglec¬ 
ting  csv  in  equation  (7U)  since  values  were  not  available,  computations  of 
Vu  were  also  made  with  the  available  results  reported  in  Ref.  (1).  All 
the  results  are  presented  in  Fig.  A36.  For  plate  D  the  predicted  Vu  are 
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in  general  larger  than  the  experimental  values  while  for  plate  B 
about  2Q)Jo  larger.  Taking  into  account  the  extreme  sensitivity  of 
equation  (7U)  to  the  values  of  (o^2  -  c^2) 2  and  csv  the  agreement  between 
predicted  and  experimental  values  of  Vu  may  be  considered  fairly  satis¬ 
factory. 


From  an  engineering  point  of  view  a  rough  estimate  of  the  ampli¬ 
tudes  in  the  velocity  range  II  can  then  be  obtained  by  approximating  the 
co0,  V  relationship  with  a  straight  line  within  this  range  as  is  indicated 
in  Fig.  15,  and  by  employing  equation  (37)  and  the  experimentally  ob¬ 
tained  values  of 


Fig.  15  Form  of  Frequency -Velocity  Relationship  in  Range  II 


The  amplitude  at  the  natural  frequency  of  the  system  is  given 
according  to  equation  (37)  by: 


a 

o 


Cln  |  H(2b)  V 

C  CO 

sv  n 
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(75) 


and  it  is  determined  uniquely  by  V,  csv,  cc^.  For  spring  II  the  velocity 
at  the  corresponding  is  larger  than  for  spring  I,  nevertheless  the 
maximum  amplitudes  for  spring  II  are  in  all  cases  smaller  than  for 
spring  I.  This  indicates  that  the  value  for  csv  may  increase  with  in¬ 
creasing  velocity  and  amplitude  of  motion.  While  csv  is  small  in  the 
present  test,  compared  with  the  Theodorsen  damping  term,  it  is  of  unique 
importance  in  the  immediate  vicinity  of  resonance  where  it  determines  the 
maximum  amplitudes. 

The  presence  of  the  third  harmonic  in  the  vibration  signal  of 
Ranges  II  and  III,  its  dependence  upon  flow  velocity  and  the  regions  of 
variable  vibrational  amplitude  are  additional  evidence  that  the  pheno¬ 
menon  under  study  is  actually  non-linear,  at  least  when  the  amplitudes 
become  large.  The  linear  model  employed  in  the  present  analysis  cannot 
allow  these  effects.  It  is  not  surprising  therefore  that  attempts  to 
use  this  linear  model  to  define  the  amplitude  and  frequency  as  unique 
functions  of  velocity  and  system  constants  are  only  moderately  success¬ 
ful  in  Range  II. 

Other  effects  such  as  non-linearities  of  the  spring  may  also  be 
present  for  large  amplitudes  and  in  spite  of  their  small  significance 
in  the  shaping  of  the  whole  response  curve  can  greatly  affect  the  maxi¬ 
mum  amplitudes  at  resonance. 

11.  Range  III  -  The  Velocity  Range  of  Small  Varying  Amplitude 
Vibrations 


Thus  far  the  discussion  and  analysis  of  the  results  have 
been  restricted  to  the  velocity  ranges  I  and  II.  The  trend  of  experi¬ 
mental  data  indicates  that  in  range  III  the  frequency  tends  to  acquire 
the  values  predicted  by  the  initial  constant  Strouhal  number.  This  has 
been  reported  by  Ippen  et  al.  (1)  but  within  the  range  of  measurable 
vibrations  in  the  present  tests  this  does  not  occur.  The  Strouhal  num¬ 
ber  increases  above  the  minimum  value  obtained  at  tc^,  V  but  still  i?  much 
smaller  than  the  initial  constant  value.  It  is  of  interest  to  study 
in  some  detail  the  moment  coefficient  diagram  for  an  understanding  of 
this  behavior.  (See  Fig.  16). 

It  is  seen  that  point  d  does  not  correspond  to  point  b.  Thus,  range  III 
is  a  transitional  range  corresponding  to  the  interval  db  where  the  fre¬ 
quency  is  governed  by  the  dynamics  of  the  system  as  in  range  II,  rather 
than  by  the  velocity  according  to  constant  Strouhal  relationship. 

Once  again,  however,  for  engineering  purposes,  it  may  be  assumed 
that  when  Vu  is  reached,  the  frequency  will  jump  immediately  to  the 
value  given  by  the  constant  initial  Strouhal  number  and  the  velocity  Vu. 
This  is  indicated  by  the  solid  vertical  lines  at  Vu  in  Figs.  A7  and  A9. 

The  Strouhal  number  then  retains  this  constant  value  until,  with  in¬ 
creasing  velocity,  the  second  harmonic  is  approached  and  the  whole  process 
repeats  itself. 


12.  The  Influence  of  Ambient  Velocity  on  the  Viscous  Damping 

In  all  of  the  previous  development  the  damping  coefficient 
in  still  water  was  considered  to  represent  the  structural  plus  viscous 
damping  of  the  system.  It  was  assumed  that  the  viscous  damping  is  in¬ 
dependent  of  the  velocity  V  and  the  motion  of  the  plate.  In  order  to 
check  the  influence  of  the  velocity  V  on  the  viscous  damping,  the 
following  tests  were  performed  with  plate  B  -  spring  II.  It  was  possible 
to  perform  free  oscillation  tests  with  the  water  running  in  the  water 
tunnel  at  the  atmospheric  pressure.  The  natural  frequency  and  damping 
coefficient  in  still  water  and  at  two  velocities  smaller  than  the  velo¬ 
city  at  which  vibrational  excitation  starts,  were  determined  with  free 
oscillation  tests.  The  results  obtained  are  shown  below. 


v  fps 

fn  cps 

Cgv/ ci 

0 

80.0 

0.9Q?o 

3.65 

80.3 

0.90?o 

h.92 

81.2 

1.03?o 

The  few  values  of  damping  coefficient  obtained  show  no  systematic  de¬ 
pendence  on  the  velocity  V  insofar  as  it  could  be  determined  experi¬ 
mentally.  The  indicative  small  increase  in  natural  frequency  with 


increasing  velocity  is  much  larger  than  that  indicated  by  inclusion  of 
the  Theodorsen  terms  and  is  within  the  experimental  error,  thus  no  de¬ 
finite  conclusion  can  be  drawn. 

An  effort  was  made  to  perform  free  oscillation  tests  at  veloci¬ 
ties  well  above  the  range  of  vibrational  excitation  of  plate  B.  No  con¬ 
sistent  results  were  obtained,  the  motion  being  irregular  and  of  high 
frequency.  From  these  limited  tests  it  appears  that  for  small  amplitudes 
and  velocities,  the  natural  frequency  and  damping  coefficient  in  still 
water  are  the  correct  values . 

13.  The  Self -Excited  Nature  of  the  Vibrations 


Analysis  of  the  experimental  results  has  indicated  a  strong 
interaction  of  forcing  moment  and  plate  motion.  Increasing  moments  are 
associated  with  increasing  amplitudes}  thus  the  vibrations  of  the  plates 
can  be  characterized  as  self-excited.  Several  phenomena  can  contribute 
to  such  strong  dependence  of  the  moment  coefficient  on  the  motion  of  the 
plate . 


Discharge  of  Vorticity 

For  the  stationary  plate  the  rate  of  discharge  of  vorticity 
in  the  wake  from  each  side  of  the  plate  is  approximately: 


> 


X-^  being  a  proportionality  constant.  For  the  oscillating  plate  with 
a  “  a0  sin  cot,  the  fluid  particles  in  the  boundary  layer  acquire  a  trans¬ 
verse  velocity,  v  0  x 'a0co  cos  cot,  in  which  x1  is  the  distance  of  the  fluid 
particle  in  the  boundary  layer  from  the  axis  of  oscillations.  Thus  the 
fluid  particles  acquire  an  additional  periodic  vorticity  component 
(3v/3x ')  **  a0 co  coscot.  The  rate  of  discharge  of  this  vorticity  in  the 
wake  is 


(  6' 

0V 

^ — r  u  dy  =  X0a  co  V  cos  cot, 

3x'  J  2  0  ’ 

-'o 

X2  being  a  proportionality  constant.  The  total  rate  at  which  vorticity 
is  discharged  in  the  wake  is  then  X^V2  +  \2d0  co  V  coscot.  A  part  of  this 
vorticity  is  diffused  and  annihilated  with  vorticity  of  opposite  sign 
and  the  rest  contributes  to  the  formation  of  discrete  vortices.  The 
presence  of  the  term  X2dQa>V  cos  cot  may  thus  be  the  cause  of  vortices  of 
larger  strength  as  a  result  of  the  plate  motion. 
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It  is  interesting  to  note  that  the  intercept  of  the  experimental 
moment  coefficient  for  a0  *  0  affords  an  estimate  of  the  lift-prottucing 
circulation  about  the  plate  in  the  absence  of  self -excitation.  From 
Eq.  (53) 


Th-UCijbV  (76) 

Assuming  that  the  vortex  is  "shed"  when  l'b(t)  *•  and  that  the  initial 

vortex  strength,  k0,  is  commensurate  with  this  circulation: 


k 


0 


rh  7.  k  ciT  b  v 

°max  a  __b_  _  -‘•I 

2n  v/Sn’  V?" 


(77) 


In  Table  III  of  Ref.  (1*),  values  of  k^^^  (i.e.,  k  at  1.5  plate 
thicknesses  downstream  of  the  trailing  edge)  are  given  for  a  velocity  of 


9.30  fps.  as  determined  from  wake  measurements  for 
vibration  prevented.  A  comparison  of  these  values 
edge)  as  given  by  Eq.  (77)  is  of  interest: 

these  same  plates  with 
with  kQ  (k  at  trailing 

ESTIMATED  VORTEX  STRENGTHS  AT  V  =  9.30  fps. 

PLATE 

ch 

k0  (trailing  edge),  ft2/sec. 
(From  Eq.  77) 

ki.5.(xA  *=i.5),ft2/s 
(From  Wake  Meas.  Ref. ■ 

A 

- 

- 

0.00127 

B 

0.0025 

0.00175 

0.0021*5 

C 

- 

- 

0.00117 

D 

0.0025 

0.00175 

0.00251 

D 

0.00125 

0.00088 

0.00251 

Three-Dimensional  Effects 


Theoretical  models  for  the  determination  of  vortex  induced 
forces  acting  on  a  body  are  usually  based  upon  two-dimensional  consi¬ 
derations  (von  Karman  vortex  street).  It  has  to  be  understood,  however, 
that  the  total  moment  acting  on  the  plate  is  an  integral  quantity  and 
that  three-dimensional  effects  can  greatly  affect  its  magnitude. 

Roshko  (31)  reports  on  the  three-dimensional  structure  of  the  wake  be¬ 
hind  a  cylinder.  His  measurements  indicate  that  the  stable  vortex 
street  (E  <  150)  has  a  periodic  spanwise  structure  with  a  wave  length 
parallel  to  the  cylinder  of  18  diameters  at  E  *  80.  Unfortunately  no 
available  results  exist  in  the  voluminous  wake  literature  concerning  the 
three-dimensional  structure  of  wakes  at  higher  Reynolds  number. 
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CD 


The  three-dimensional  effects  are  of  extreme  importance  as  is 
shown  tfelow. 


The  forcing  moment  per  unit  spanwise  length  of  the  plate  at  any 
spanwise  position,  z,  can  be  expressed  as 

mf  l  p(2b)2V2  e“*. 

The  coefficient,  mf  ,  as  determined  by  measurement  of  the  net  moment  on 

the  entire  plate  is°then  given  by:  „ 

f  h 


mfr 


i 


dz 


If  m^  is  periodic  such  that  H  is  an  integral  multiple  of  the 
spanwise  wavezlength  then  mf0  ■=  0. 

If  mf  is  periodic  but  H  is  not  an  integral  multiple  of  the 
spanwise  wave  length  then  mf  may  have  any  value  up  to  the  spatial 
RMS  of  mf  .  0 

For  the  stationary  plate,  the  degree  of  two-dimensionality  of 
the  wake  should  depend  strongly  on  body  Reynolds  number  and  perhaps  to 
some  degree  on  the  trailing  edge  geometry.  For  an  oscillating  plate, 
the  degree  of  two-dimensionality  may  logically  be  expected  to  increase 
with  the  amplitude  of  motion,  at  least  when  the  frequency  of  oscillation 
is  that  of  the  wake  motion  as  it  is  in  the  forced  vibrations  being  con¬ 
sidered  here.  The  observed  increase  of  moment  coefficient  with  increase 
in  the  dynamic  parameter: 

„  2bcoa 

T7-*  _  O 


may  thus  result  from  increases  in  two-dimensionality  as  well  as  from  in¬ 
creases  in  local  vortex  strength. 

Plates  B  and  D  with  extending  trailing  edges  exhibit  large  ampli¬ 
tudes  and  thus  large  values  of  V*  in  the  vicinity  of  resonance.  It  is 
thus  possible  that  the  extension  of  trailing  edges  beyond  the  separation 
points,  and  thus  the  provision  of  solid  shielding  between  the  shearing 
layers  of  opposite  sign,  greatly  affects  the  strength  of  the  shed  vor¬ 
tices  in  addition  to  the  two-dimensional  structure  of  the  wake.  For 
plates  A  and  C  the  trailing  edge  does  not  provide  solid  shielding  of  the 
shearing  layers  and  larger  annihilation  of  vorticity  may  take  place  in 
the  early  wake,  resulting  in  vortices  of  smaller  strength. 


D.  Final  Tests  -  Wake  Measurements 


1.  Experimental  Procedure 

For  several  plate-spring  systems,  several  discrete  veloci¬ 
ties  were  selected  within  the  velocity  range  of  measurable  vibrational 
excitation  and  at  each  velocity  the  amplitude  and  frequency  of  vibration 
were  measured.  Measurements  of  the  RMS  of  the  total  head  fluctuations, 
H't,  were  then  made  across  the  half -width  of  the  wake  at  station 
x  ■  l.£t  utilizing  the  thermocouple  voltmeter.  Measurements  were  taken 
at  sufficient  transverse  points  so  that  a  continuous  distribution  curve 
could  be  drawn.  Three  measurements  of  the  RMS  were  made  at  every  point 
and  the  average  of  them  considered  as  the  true  RMS  value.  The  discrete 
frequencies  of  the  spectrum  of  't2  were  picked  up  with  the  vibration 
analyzer. 


The  plate-spring  systems  and  the  number  of  velocities  for  which 
transverse  distributions  of  the  total  head  fluctuation  were  made  are  in¬ 
dicated  below. 


Plate  B  -  spring 

I 

7  Runs 

Plate  C  -  spring 

I 

3  Runs 

Plate  C  -  spring 

II 

3  Runs 

Plate  D  -  spring 

I 

10  Runs 

Plate  D  -  spring 

II 

U  Runs 

2.  Interpretation  of  Wake  Measurements 

The  response  of  the  turbulence  gage  used  for  the  fluctuating 
total  head  measurements  is  discussed  in  detail  in  Refs.  (2),  (3)  and  (I4.3) . 
It  is  shown  in  these  references  that,  neglecting  the  probe  response  to 
the  transverse  components  of  turbulence,  the  RMS  of  the  fluctuating  com¬ 
ponent  of  the  local  instantaneous  total  head  is  given  by: 


in  which 

em  =  amplified  gage  signal  input  to  the  measuring  instruments 
Ga  =  amplifier  gain 
Gp  B  preamplifier  gain 

Kq  ■  gage  sensitivity  in  millivolts  per  foot  of  water 
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Under  additional  restrictive  conditions,  the  principle  one  being 
that  Vu*2  «  u,  Ippen  (1*2)  has  shown  that  Eq.  (78)  may  be  simplified  to 


which  provides  a  means  of  studying  the  longitudinal  turbulent  velocity 
component,  u',  provided  the  mean  velocity,  u,  is  known. 

Under  the  conditions  existing  in  the  early  wake,  where  the  pre¬ 
sent  measurements  were  made,  the  applicability  of  Eq.  (79)  is  very 
questionable.  A  reduction  of  data  from  the  preliminary  tests  performed 
with  plate  C  -  spring  land  using  Eq.  (79)  indicated  local  values  of  the 
relative  turbulence  N/u'2  / u  as  high  as  1  on  the  wake  centerline  and 
about  0.1  at  a  transverse  distance,  y/t  ■  1.1.  For  this  reason,  no 

attempt  to  calculate  Vu^  from  this  reduced  relation  was  felt  justi¬ 
fied  and  mean  local  velocity  traverses  were  not  made.  The  reduced 
formula  will  be  used  only  for  a  general  qualitative  interpretation  of 
the  measurements  taken. 

The  actual  \J H1^  consists  of  two  parts:  a  random  turbulent 

background  covering  the  continuous  spectrum  of  frequencies  and  discrete 
components  due  to  the  "shedding  of  vortices".  In  these  tests  the  dis¬ 
crete  component  at  the  fundamental  vibrational  frequency  accounts 
usually  for  more  than  of  the  total  RMS  as  has  been  discussed  in 

section  V-B-3. 

For  an  ideal  vortex  street  (von  Karman  vortex  street)  with  geo¬ 
metric  configuration  independent  of  the  ambient  velocity  and  of  plate 
motion,  the  N/u  * z  at  any  particular  point  in  the  wake  will  be  pro¬ 

portional  to  the  strength  of  the  discrete  vortices. 

N/u,a  'v? vortex  strength 

Under  the  assumption  of  similar  velocity  distribution  inthe  wake  it  can 
be  deduced  that  for  a  given  x/t  in  the  early  wake  u  v/u'2  ~  V \/u 12  ~V» 
(Vortex  strength)  V •  (circulation)  ~  lift  force  ^  vortex-induced  lift 
moment.  For  the  qualitative  arguments  of  the  present  discussion  we  will 
assume  that 


u 


in  spite  of  the  inaccuracy  pointed  out  above.  Thus  it  is  concluded  that 

the  sJW?  at  any  particular  point  in  the  wake  is  very  roughly  propor¬ 
tional  to  the  lifting  moment  acting  on  the  plate. 
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To  avoid  reliance  upon  a  point  measurement  for  comparative 
studies  the  integral  of  n/h1^®  d(y/t)  over  the  width  of  the  wake 
will  be  used  as  a  better  index  of  a  two-dimensional  vortex- induced 
lifting  moment  than  is  \/h  '^a  at  a  single  location. 


3.  Correlation  of  Wake  Measurements  with  Vibrational 
Characteristics 


The  transverse  distributions  of  the  total  head  fluctuation  in 
the  wake  of  the  plates  (at  midspan  only)  are  presented  in  Figs.  A37 

through  AltO.  The  dimensionless  total  head  fluctuation  (2g  ")AZ 
has  been  plotted  versus  the  dimensionless  transverse  distance  y/t  over 
the  half  width  of  the  wake,  6.  The  quantity 


m 
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has  been  computed  for  every  test  velocity  by  numerical  integration.  The 
integral  has  been  extended  out  to  y/t  “1.1  since  it  was  found  in  the 
preliminary  tests  that  there  the  mean  local  velocity  is  almost  equal  to 
the  free  stream  velocity.  The  computed  values  of  m  as  well  as  the 
ambient  velocity,  amplitude  and  fundamental  vibrational  frequency  are 
presented  in  the  Figs.  A37  through  Ai*0. 

It  is  obvious  from  the  developments  of  the  previous  section 
that  the  dimensionless  quantity  m  should  be  roughly  proportional  to  the 
two-dimensional  forcing  moment  coefficient.  Also,  the  effect  of  in¬ 
creasing  two-dimensionality  of  the  wake  cannot  be  reflected  in  the  co¬ 
efficient  m  due  to  the  single  spanwise  location  of  the  measurements 
made.  A  study  of  the  values  of  m  reveals  the  following: 

For  plate  D  (Fig.  A3 9  )  the  value  of  m  remains  essentially  con¬ 
stant.  The  extent  of  the  influence  of  the  plate  motion  on  the  strength 
of  the  vortices  is  illustrated  best  perhaps  by  Fig.  AhO  in  which  we  can 
compare  values  of  m  for  different  a0  but  at  the  same  velocity.  For  a 
3-fold  variation  in  amplitude  (V  =  0.37  fps),  m  increases  by  a  factor  of 
only  1.6  yet  from  Fig.  k}k  the  moment  coefficient  increases  by  a  factor 
of  2.7.  We  are  thus  led  to  conclude  that  the  sensitivity  of  moment  co¬ 
efficient  to  plate  motion  is  at  least  as  dependent  upon  other  mechanisms 
such  as  increasing  wake  two-dimensionality  or  transverse  redistribution 
of  vorticity  as  it  is  upon  local  increase  in  the  strength  of  discrete 
vorticity. 


For  plate  B  (Fig.  A37),  m  undergoes  a  1.7  fold  variation.  Again 
using  the  corresponding  values  of  the  self -excitation  parameter,  V*,  the 
moment  coefficient  is  seen  to  increase  by  a  factor  of  7.2.  For  this 
plate,  the  increase  in  local  vortex  strength  is  clearly  a  minor  effect. 

For  plate  C  (Fig.  A38)  m  changes  by  a  factor  of  1.26  while  mf0 
experiences  only  a  1.16  fold  variation.  Due  to  the  qualitative  nature 
of  the  total  head  interpretations,  speculations  concerning  the  sign  of 
this  difference  between  m  and  mfQ  are  unwarranted  but  it  appears  as 
though  increasing  vortex  strength  may  account  for  the  entire  increase 
in  moment  coefficient. 

Although  the  above  comparisons  are  only  qualitative,  the  nature 
of  the  major  inaccuracies  in  interpretation  of  the  probe  response  are 
such  as  to  give  turbulent  velocities  which  are  too  large.  For  this 
reason,  the  indication  that  at  least  certain  of  the  plates  become  self- 
excited  through  increasing  wake  two-dimensionality  bears  emphasis  by 
a  tabulation  of  the  factors  just  discussed: 


Importance  of  Increasing  Local  Vortex  Strength 
in  the  Self-Eixcitation  Mechanism 


Plate 

Maximum  Factor  of 
Variation  in  "m" 

Corresponding  Factor  of 
Variation  in  mf0 

Apparent  Major 

Self -Excitation  Mech. 

B 

1.7 

7.2 

Two-dimensionality 

C 

1.6 

2.7 

Two-dimensionality 
and  local  vortex 
strength 

D 

1.26 

1.16 

Local  vortex  strength 

It  should  also  be  noted  that  for  all  transverse  positions  dis¬ 
crete  frequencies  were  present  in  the  spectrum  of  _  ^he  fundamen¬ 

tal  frequency  was  easily  detected  as  were  the  second  and  third  harmonics. 

The  fundamental  discrete  frequency  in  the  spectrum  of  \/h  't2  was  exact¬ 
ly  the  same  as  the  fundamental  frequency  of  the  vibrating  plate  in  all 
cases. 


1:.  The  Periodic  Flow  Field  in  the  Wake 


The  electrical  signal  of  H'^  was  observed  on  the  screen  of 
the  oscilloscope.  Near  the  centerline  of  the  wake  the  signal  possessed 
periodicity  but  it  was  obscured  by  the  heavy  turbulent  background.  The 
periodicity  was  more  clear  as  the  probe  was  moved  away  from  the  center- 
line.  From  the  investigation  of  the  wake  structure  for  stationary 
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plates  (2) (3) ,  It  is  known  that  the  probe,  at  large  values  of  y/t, 
responds  to  the  free  turbulence  level.  No  discrete  periodicities  are 
present  at  this  position  with  the  stationary  plates  and  the  spectrum 

of  n/h*^  is  continuous.  In  the  present  case  of  the  vibrating  plate 
the  signal  of  H’^  at  the  points  y/t  -  1.82,  2.62,  5.1*2  was  periodic 
and  almost  sinusoidal  with  a  very  small  turbulent  background  superim¬ 
posed.  This  was  observed  for  plates  B  and  D  which  give  large  vibra¬ 
tional  amplitudes .  The  RMS  level  of  H'^  at  these  points  is  consider¬ 
ably  higher  than  the  corresponding  level  of  thefree  turbulence.  The 
free  turbulence  level  in  the  water  tunnel  is  *  3?o  or  l*/o 

while  the  RMS  of  the  H1^  at  y/t  ■  5.1*2  gives  an  indicated  value  of 
*  25^o  for  plate  B  at  resonance. 

A  special  test  was  carried  out  with  plate  B  -  spring  I  when 
this  observation  was  made.  The  total  head  fluctuation  gage  was 
placed  at  station  x/t  “1.5,  y/t  =  5.1*2 .  The  velocity  in  the  water 
tunnel  was  increased  by  small  increments  and  the  RMS  of  H1^  was 
measured  by  the  A.C.  Voltmeter.  The  test  covered  the  whole  range 
of  velocities  where  vibrational  excitation  occurs.  The  vibrational 
amplitude  versus  velocity  relationship  having  been  established,  a 
direct  comparison  could  be  made  between  w  and  the  vibrational 

amplitude;  It  was  discovered  that  for  large  amplitudes  the  n/hmF 
at  y/t  “  5.U2  is  proportional  to  the  vibrational  amplitude.  The  same 
behavior  is  exhibited  with  plate  D  for  large  amplitudes. 

Table  VIII  presents  the  values  of  at  y/t  ■  5.1*2,  the 

corresponding  velocity,  V,  the  RMS  of  vibrational  amplitude,  oQ,  and 

the  ratio  o.0/sj H for  plates  B  and  D. 

For  plate  B  the  values  of  \/^*t2  at  the  transverse  positions 
y/t  *  1.82,  2.62,  5.1*2  remain  almost  constant  with  higher  values  some- 
tlmes  at  y/t  ”  5.1*2  as  it  is  shown  in  Fig.  A3 7.  For  plate  D  the 

x/S't5  decreases  slightly  from  y/t  ■  1.82  to  y/t  -  5.1*2.  For  plate  C 

the  value  of  v/H'^a  at  y/t  *  1.82,  2.62,  5.1*2  is  larger  than  the  free 
turbulence  level  but  the  periodicity  is  obscured  by  the  turbulent  back¬ 
ground.  The  points  y/t  ■  1.82,  2.62,  5.1*2  at  the  station  x/t  ■  1.5  are 
outside  of  the  conventionally  described  wake  width. 

From  the  complete  expression  for  the  total  head  fluctuation  it 
is  observed  that  two  terms  can  make  a  significant  contribution,  u'  and 
p’.  It  is  thus  likely  that  a  fluctuating  pressure  field  is  created  by 
the  vibrating  plate  with  p'  a  ■  o0  cos  co0t.  In  this  case 
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TABLE  VIII 


Free  Stream  Total  Head  Fluctuations 
and  Vibrational  Amplitudes 

Plate  B  -  Spring  I 


10^  rad 

at  y/t  ■ 

ft.  water 

5.1*2  0 

10J  rad 

V 

fps 

2  c . 

*  ft.w. 

5.58 

1.20 

0.074 

16.2 

5.67 

2.50 

0.104 

24 

6JI4.O 

8.70 

0.344 

25.3 

6:70 

9.25 

0.382 

24.2 

6.82 

9.38 

O.388 

.  24.2 

Measurements  of 

7.13 

7.50 

9.05 

8.60 

0.374 

0.333 

24.2 

25.8 

x/h1*.2  made  with 
A.C.  Balantlne 

7.88 

8.00 

0.317 

25.2 

Voltmeter. 

8.27 

7.35 

0.294 

25.0 

8.80 

0.90 

0.118 

7.6 

8.16 

0.35 

0.117 

3.0 

5.7? 

3.61 

0.131 

27.5 

6.00 

6.43 

0.245 

26.2 

Measurements  of 

6.47 

7.17 

7.50 

9.15 

8.91 

7.98 

0.301 

0.358 

O.320 

30.4 

24.9 

24.9 

Vh  't2  made  with 
Thermocouple 

V  oltmeter # 

8.60 

1.08  2.35 

0.216 

10.9 

Plate 

D  -  Spring  I 

6.37 

0.86 

0.084 

10.2 

8.00 

3.15 

0.173 

18.2 

9.29 

4.87 

0.245 

19.85 

10.80 

7.02 

0.353 

19.9 

11.75 

7.92 

0.353 

22.4 

Measurements  of 

12.40 

13.65 

14.20 

15.05 

8.64 

8.89 

8.74 

8.06 

0.431 

0.465 

0.490 

0.490 

20.0 

19.1 

17.8 

16.45 

made  with 
the  Thermocouple 
Voltmeter. 

15.50 

6.83 

0.490 

13.9 

Plate 

D  -  Spring  II 

7.97 

2.17 

0.142 

15.3 

10.8 

5.39 

0.314 

17.1 

12.45 

6.74 

0.441 

15.3 
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As  long  as  p^/y2  »  uz  uP^/g2,  the  \J H'^2  at  y/t  -  5^42  is  propor¬ 
tional  to  Op.  The  importance  of  the  term  u2  u* Vg2  increases  with  in¬ 
creasing  velocity.  The  above  arguments  are  greatly  supported  by  the 
experimental  results. 

A  careful  study  of  Table  VIII  reveals  the  following.  For  plate  B 
in  the  small  range  of  relatively  low  velocities  where  large  amplitudes 

are  obtained,  the  proportionality  of  to  the  vibrational  amplitude 

is  perfect,  because  the  term  u2  u,2/g2  is  of  minor  importance  compared 
with  p '  2/y2 .  To  the  left  and  right  of  resonance,  with  small  amplitudes 

and  thus  increasing  importance  of  the  term  u2  u,2/g2,  the  proportionality 
fails.  For  plate  D,  which  exhibits  large  amplitudes  at  relatively  high 
velocities,  the  importance  of  the  term  u2  ul2/g2  is  greater  than  for  plate 
B  and  the  proportionality  of  \Jyl  1  ^  to  the  vibrational  amplitude  holds 
only  near  resonance. 

To  see  if  the  periodic  pressure  field  created  by  the  plate 
vibration  is  imposed  upon  the  whole  flow  field,  the  total  head  fluctua¬ 
tion  gage  was  placed  at  station  x/t  -  10  and  y/t  -  5.U2.  Ttie 

at  this  position  was  found  to  be  only  2/3  of  the  value  obtained  at 
x/t  -1.5  and  y/t  -  5«it2.  The  signal  of  H't  exhibited  random  charac¬ 
teristics  with  an  almost  complete  absence  of  discrete  periodicity.  This 
observation  indicates  that  the  free-stream  periodicities  are  likely  the 
result  of  fluctuating  local  pressures  caused  by,  and  related  in  magni¬ 
tude  to,  the  local  transverse  velocity  fluctuations.  These  latter  fluc¬ 
tuations,  near  the  plate,  will  be  directly  related  to  the  magnitude  of 
aQ,  while  with  distance  downstream  in  the  wake  they  are  subject  to  vis¬ 
cous  decay  and  degeneration  into  random  turbulence. 

It  was  thought  of  interest  to  investigate  the  phase  difference 
between  the  motion  of  the  plate  and  the  periodic  flow  disturbance  out¬ 
side  the  wake.  For  this  purpose  a  special  test  was  performed  with  plate 
D  -  spring  II.  The  total  head  fluctuation  gage  was  placed  at  x/t  -  1.5 
and  y/t  -  5»i*2.  The  signals  H1^  and  a  were  connected  to  the  two  channels 
of  the  oscilloscope  and  were  observed  on  the  screen.  For  all  tested 
velocities  the  signals  were  in  phase  as  far  as  it  could  be  determined  ex¬ 
perimentally.  Fig.  Aitl  presents  pictures  of  the  two  signals  under  two 
of  the  test  conditions.  The  observed  phasing  is  what  would  be  expected 
from  the  dynamic  pressure  changes  accompanying  an  oscillating  velocity 
field.  That  is,  the  pressure  should  be  ?0°  out  of  phase  with  the  fluid 
velocity  and,  close  to  the  plate,  the  fluid  velocity  should  be  90°  out  of 
phase  with  the  plate  motion  which  initiates  it. 
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71.  SUMMARY  AUD  CONCLUSIONS 


The  vibrational  behavior  and  wake  structure  of  four  flat  plates 
(chord  -  2",  thickness  ■  l/li")  with  different  trailing  edge  geometry, 
mounted  at  zero  mean  angle  of  attack  in  a  uniform  flow  field  and  allowing 
torsional  oscillations  about  an  axis  containing  the  leading  edge,  has 
been  investigated  in  a  closed- jet  water  tunnel. 

Measurements  of  the  amplitude  and  frequency  of  the  vibration  and 
of  the  fluctuating  total  head  in  the  wake  were  made  over  the  range  of 
ambient  velocity  producing  measurable  vibrational  excitation. 

An  equation  for  the  amplitude  response  of  the  plate-spring  system 
is  developed  incorporating  the  hydrodynamic  loads  predicted  by  the  linea¬ 
rized  potential  theory  of  Theodorsen  and  an  unknown  flow- induced  forcing 
moment. 


Using  the  experimental  data  and  the  above  theoretical  model  the 
forcing  moment  coefficients  were  evaluated  and  correlated  with  the  dyna¬ 
mic  characteristics  of  the  flow-plate  system. 

The  following  conclusions  are  drawn: 

1.  Similarity  of  normalized  spectral  distributions  of  total  head 
fluctuations  in  the  early  wake  of  the  same  plate  when  stationary  and  when 
vibrating  indicates  that  periodic  vortex  shedding  provides  the  forcing 
mechanism  in  the  flow- induced  plate  vibration  studied. 

2.  The  fundamental  vibrational  frequency  coincides  exactly,  in 
all  cases,  with  the  fundamental  discrete  frequency  in  the  total  head 
fluctuation  spectrum. 

3.  Large  amplitude  flow-induced  vibrations  are  characterized  by 
three  distinct  ranges  of  vibrational  behavior: 

Range  I:  "simple  harmonic"  motion  with  slowly  varying  small 
amplitude  where  a  constant  body-thickness  Strouhal 
No.  (St  =  f0t/7)  prevails.  This  range  ends  at  a 
frequency  less  than  the  resonant  frequency. 

Range  II:  periodic  motion  with  large  stable  amplitude  con¬ 
taining  a  third  harmonic  which  grows  with  ambient 
velocity.  Strouhal  number  (based  on  the  fundamen¬ 
tal  frequency)  continually  decreases  with  ambient 
velocity  in  this  range  which  ends  at  a  frequency 
higher  than  the  resonant  frequency. 
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Range  Ills  periodic  motion  of  small,  unstable  amplitude. 

Strouhal  number  continually  increases  with 
ambient  velocity. 

U.  The  Range  I  Strouhal  number  appears  to  be  the  same  as  the 
steady  body  Strouhal  number. 

!>.  Vortex-induced  forcing  moments  and  plate  motion  are  related 
uniquely  for  all  plates  tested  through  the  self-excitation  parameter, 

V*  "  2bto0a0/V .  This  parameter  is  the  ratio  of  the  RMS  transverse  velo¬ 
city  of  the  trailing  edge  to  the  mean  ambient  flow  velocity.  It  is 
shown  to  be  proportional  to  the  ratio  of  separation  stream  line  velocity 
to  free  stream  velocity  and  hence  is  also  proportional  to  wake  pressure 
coefficient. 


The  form  of  this  relationship  is  a  piecewise  linear  variation  in 
the  forcing  moment  coefficient,  mf0,  with  the  self -excitation  parameter, 
V*.  Corresponding  to  the  above  Ranges: 


Range  I: 

mf0  “  Ciz  +  C2-,-  V* 

(28) 

Range  II: 

% o  "  Cln  +  C2n  v* 

(29) 

6.  The  linear  dependence  of  forcing  moment  of  the  self -excita¬ 
tion  parameter  yields  the  effect  of  a  decreasing  damping  with  increasing 
velocity.  It  leads  to  an  eigenvalue  problem  defining  the  critical  fre¬ 
quencies  bounding  Range  II  in  terms  of  the  system  constants,  the  "moment 
slope",  Cj,  and  the  Range  I  Strouhal  number,  St,  i.e. 


o  2  ®  t  2 

®o  }  +(csv  +  C(2^)} 


o  ' «-  ®  "t 

“o  1  "  CI  (2nS7)a>o 


0 


(U3) 


7.  The  ambient  velocity,  Vj,,  corresponding  to  the  lower  critical 
frequency,  cc^,  (and  the  beginning  of  large  amplitude  motion)  can  be  pre¬ 
dicted  precisely  through  use  of  the  smaller  root  of  Eq.  (1*3)  and  the 
Range  I  Strouhal  number,  S^. 

ox  t 

VL  =  2^  (6l) 


8.  The  ambient  velocity,  Vu,  corresponding  to  the  upper  critical 
frequency,  toy  (and  the  end  of  large  amplitude  motion)  can  be  predicted 
accurately  enough  for  engineering  purposes  (+5  to  20/o)  by  assuming  the 

reduced  amplitude ,  aQ/V 


to  be  the  same  at  and  a 
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This  yields: 


.i.i^V-02  , 

5  ‘5  (A-  c  Z~  ^ 
sv  u  u 


where 


A  -  IVV-V) 


+  (c  + 
v  sv 


CV, 


2  2 
)  “L 


% 

^  -  CTL“L 


(71*) 


(70) 


9.  Controlled  variation  of  the  system  parameters  should  allow 
design  of  a  system  in  which  Eq.  (1*3)  has  only  one  root,  i.e.  where 
Range  II  is  absent  and  thus  large  flow-induced  vibrations  are  not  in¬ 
curred.  This  can  perhaps  be  most  practically  accomplished  by  increasing 
the  combined  structural  and  viscous  damping  coefficient,  cgv,  until 

<50) 

For  a  fixed  value  cf  the  damping,  csv,  the  critical  frequency 
range,  ^ can  be  reduced  through  an  increase  in  the  Strouhal  number, 

st‘ 


10.  Wake  measurements  suggest  that  the  large  increase  in  moment 
coefficient,  rrif  ,  with  increasing  self-excitation  parameter,  V*,  may  be 
due  to  either  l8creasing  two-dimensionality  of  the  wake  or  to  increasing 
local  vortex  strength  depending  upon  the  trailing  edge  geometry.  For  the 
plates  tested  the  former  phenomenon  seems  to  govern  when  the  trailing 
edge  separation  points  are  variable  and  the  latter  when  separation  points 
are  fixed. 


VII.  APPENDICES 


APPENDIX  A:  Experimental  Results 

APPENDIX  B:  Spring  Calibration  and  Free  Oscillation  Tests 

APPENDIX  C:  References 
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Dimensionless  Transverse  Distance 


PLATE  C  E 

HU  Spring  1 

Station  -y— 

*  1.5 

o  Vibrating 

Plate 

•  Stationary 

Plate 

1.0  1.0  1.0  1.0  1.0 


Dimensionless  Velocity  -g- 


Tig.  A1  Mean  Velocity  Distributions  in  the  Wake 


(Note:  Units  of  velocity,  V,  are  fps) 


1  I  •  I  I  >  I  '  !  '  I  I  '  *  '  ' 


Fig.  JULO  Vibrational  Signal 
Plate  B  -  Spring  I 
Range  II 

Stable  Vibrations 
V  -  6.li7  fps,, 
o"  -  9.15  10  3  rad., 
f°  -  71.2  cpa. 


A 


Fig.  112  Vibrational  Signal 
Plate  B  -  Spring  I 
Point  of  Interaittancy 
between  Range  II, III 
V  -  8.60  fps,- 

-  2.35-1.08  10"  3  rad., 
f,-  76.3  eps- 


Fig.  All  Vibrational  Signal 
Plate  B  -  Spring  I 
Range  II 

Stable  Vibrations 
V_-  7.17  fps, 

Oq  ••  8.91  10“3  rad., 
fD  »  72.6  cps. 


Fig.  A13  Vibrational  Signal 
Plate  B  -  Spring  I 
?oint  of  Intemdttancy 
between  Range  II, III 
V  ■  8.60  fps, 

oS  -  2.S5-1.08  10‘3  rad., 
t9m  76.3  cps. 
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Fig.  Allt  Vibrational  Signal 

Fig.  A15  Vibrational  Signal 

Plate  D  -  Spring  I 

Plate  D  -  Spring  I 

Range  I 

Range  II 

Unstable  Vibrations 

Stable  Vibrations 

V  -  6.37  fpe. 

SJ  -  0.86  10-3  P«d., 

V  -  9.29  fps, 

do  ■  It. 86  10  -  rad., 

f0  »■  6b. It  cps. 

f0  *  69.7  cps. 

Fig.  A17  Vibrational  Signal 
Plate  D  -  Spring  I 
Range  II 

Stable  Vibrations 
V  ■  lit. 20  fps. 

■  8.7b  10-3  rad., 
fe  "  75.5  ope. 
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Fig.  A16  Vibrational  Signal 
Plate  D  -  Spring  I 
Range  II 

Stable  Vibrations 
V  -  9.29  fps, 

■  h.86  10’ 3  rad., 
f0  ■  69.7  cps. 


pic.  A18  Vibrational  Signal 
Plate  C  -  Spring  II 
Unstable  Vibrations  of 
varying  amplitude 
V  -  6.I46  fps. 

35  -  1.52  10‘3  rad., 
fD  -  7h.5  cps. 


Fig.  A19  Vibrational  Signal 
Plate  C  -  Spring  II 
Unstable  Vibrations  of 
varying  amplitude 
V  -  7-0?  fps., 
oS  •  1.75  10"3  rad., 
f0  -  77.0  cps. 
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4  5  6  7  8  9  10  II 

Ambient  Velocity  V  fps 


Pig.  A23  Reproducibility  of  Results,  Plate  B  -  Spring  II 
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RMS  of  Vibrational  Amplitude  10  rad. 
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Pig.  A26  Response  Curve,  Plate  C  Pig-  A27  Response  Curve,  Plate  D 


Fig.  A 29  Reduced  Response  Curre,  Piste  B 


Fig.  A 30  Reduced  Response  Curve,  Piste  C  Fig.  A31  Reduced  Response  Curve,  Piste  D 
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Fig.  A32  The  Moment  Coefficient 


Fig.  A33  The  Moment  Coefficient,  Expanded  Scale 
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Measured  Critical  Frequencies  ,  fL  ,  ,  cps 


"’Self  Excitation.  Parameter  \A  |qb 


Fig.  A3li  The  Moment  Coefficient,  Hysteresis  Effect 


Predicted  Critical  Frequencies,  f L  ,  f u ,  cps 


Fig.  A3<  Prediction  of  Critical  Frequencies 
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Pig.  Al6  Prediction  of  Critical  Velocities 
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Fig.  A38  Transverse  Wake  Profiles  -  Fluctuating  Total  Head,  Plate  C 
(Notes  In  above  figures,  V  in  fps,  a0  in  10" ^  rad.,  f0  in  cps) 
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Plat*  D-Spring  II 
V  -  13. SO  fps  a0  -  7.U1 10  3  rad. 

fo  "  cps 
Total  H*ad  Fluctuation 
Qag*  at  ^  -  5.U2 


Plata  D-Spring  n 

V  •  9.60  fps  a  »  U.UU10"3  pad. 
0 

f  ■  73.0  cps 

Total  H*ad  Fluctuation 
Qag*  at  ^  -  5.42 


Pig.  Aul  Vibrational  Signal  t*.  Total  Head  Fluctuation  Signal 
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APPENDIX  B 


Spring  Calibration  and  Free  Oscillation  Tests 
1.  Spring  Calibration 

It  can  be  shown  theoretically  (U3)  that  the  unbalanced  voltage, 

AE,  (in  millivolts)  across  a  U  active  arm  strain  gage  bridge  arrayed 
for  torsion  measurement,  is  related  to  the  excitation  voltage,  E,  (in 
volts)  and  the  angular  deflection,  a,  (in  radians)  by  the  expression: 

AE  ■*  G  E  a 
s 

in  which  Gg  is  the  torsion  spring  calibration  constant  in  millivolts  per 
radian  per  volt  of  excitation. 

0S  was  determined  experimentally  by  an  in-place  calibration.  To 
accomplish  this,  the  set-screws  which  fasten  the  plate  shaft  to  the  spring 
were  replaced  by  two  long  stainless  steel  rods  threaded  at  one  end  to  in¬ 
sure  the  shaft-to-spring  rigidity.  Applying  a  force  to  one  rod  caused 
the  other  to  deflect.  This  deflection  was  measured  to  the  nearest  0.001 
inch  with  a  micrometer  at  a  distance  of  10  inches  from  the  spring  axis. 

The  bridge  output  was  amplified  500  or  1000  times  and  measured  with  a 
D.C.  voltmeter. 

Each  spring  was  tested  in  both  directions  to  check  its  symme¬ 
trical  behavior  and  with  increasing  and  decreasing  deflection  to  check 
any  hysteresis.  Each  spring  was  calibrated  in  place  at  the  beginning 
and  at  the  end  of  vibrational  tests  with  every  plate-spring  system  in 
order  to  check  possible  spring  fatigue.  Fig.  B1  shows  one  of  the  cali¬ 
bration  graphs  obtained. 

The  spring  sensitivity  Gs  for  each  plate-spring  system  is  in¬ 
dicated  below. 


Exploratory  Tests 

Plate  C  -  Spring  I  Gs  =  16.9  ^-volfol  Batt. 

Final  Tests 


Plate 

B 

D 

C 

A 


Gs  rad-volt  of  Battery 


Spring  I 
17.2 
17.3k 
16.78 
16.78 


Spring  II 
18.13 
18.27 
18.26 
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Sample  of  Spring  Calibration 


Plate  D  -  Spring  II 
Beginning  of  Testa 


fn  air  "  153  CP8* 


p-  -  1.01*2% 


. in  Mil 

Ml.,,.'*'# 


Plate  B  -  Spring  I 
Beginning  of  Tests 

fR  •  75.8  cps. 

^-1.05O/o 
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B2  Free  Oscillation  Test  in  Air  B3  Free  Oscillation  Test  in  Still  Water 
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Simultaneous  measurement  of  the  force  applied  to  produce  the  de¬ 
flections  allows  determination  of  the  structural  spring  constant,  Ks. 

2.  Free  Oscillation  Tests 

In  order  to  determine  the  plate-spring  system  constants,  free  os¬ 
cillation  tests  were  performed  in  air  (water  tunnel  empty)  and  in  still 
water : 


Free  oscillation  tests  in  air: 

Assuming  linear  damping,  the  equation  of  motion  for  the  plate-spring 
system  in  air  is: 

Ia+ca+Ka“0 
s  s  s 

Under  the  initial  conditions: 


a  =  a  and  a  1  0  at  t  *  0, 
o  * 

the  solution  of  the  above  equation  yields 


/  K  c  2 

co  in  air  *  - 

n  ' 1  h  Uit2 


or  for 


1*1+ 


« 


r  : 

s 


“n 111  a*T  ~\Jrs 

Free  oscillation  tests  in  still  water: 

Again  assuming  linear  dairying,  the  equation  of  motion  for  the  plate-spring 
system  in  still  water  is: 

(I  +  I  )  a  +(c  +c)a+Ka“0 
v  s  p  '  s  v'  s 


or 


I.  o+c  a  +  K  a  *  0 
t  sv  s 


Under  the  initial  conditions: 


o  *  a  and  a  ■  0  at  t  ■  0, 
o 
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the  solution  yields 


co  in  still  water 
n 


or  for 


Ks 


li  “ 

/  s  sv 


co  in  still  water 
n 


For  the  tests  in  air  the  plate  was  very  easily  deflected  and  suddenly 
released  by  striking  it.  The  decaying  oscillation  was  photographed  by 
using  a  polaroid  camera  fitting  on  the  screen  of  the  oscilloscope.  For 
the  tests  in  still  water  it  was  impossible  to  strike  the  plate;  instead 
the  spring  itself  was  given  a  deflection  and  then  suddenly  released. 

Figures  B2,  B3  are  samples  of  the  free  damped  oscillations  ob¬ 
tained. 


Knowing  the  sweep  rate  of  the  oscilloscope  signal  the  frequency 
of  oscillation  was  easily  obtained  from  the  pictures.  Several  pictures 
were  taken  for  each  plate  spring  system  and  the  average  value  of  fre¬ 
quency  fn  in  air  and  fn  in  still  water  was  considered  the  true  value. 

The  experimental  error  in  determining  the  frequency  from  the  pictures  is 
about  1.5/0. 

The  damping  coefficients  cs  and  cgv  were  also  obtained  from  the 
free  oscillation  pictures  by  the  conventional  method  of  determining  the 
average  logarithmic  decrement  for  successive  oscillations.  For  each 
plate-spring  system,  the  coefficients  were  determined  from  several  pic¬ 
tures  and  averaged.  Although  the  assumption  of  linear  damping  proved 
satisfactory,  the  individual  coefficients  showed  a  scatter  of  up  to 
lO/o  of  the  average  value. 

During  the  course  of  experiments  it  was  discovered  that  the 
frequency  recorded  by  the  vibration  analyzer  and  the  oscilloscope  were 
not  exactly  the  same.  Therefore  a  calibration  test  of  the  oscilloscope 
versus  the  vibration  analyzer  was  performed.  A  sinusoidal  electrical 
signal  from  a  signal  generator  was  fed  to  the  vibration  analyzer  and  the 
oscilloscope  simultaneously.  From  pictures  taken  the  frequency  given  by 
the  oscilloscope  was  determined  and  plotted  versus  the  frequency  obtained 
by  the  vibration  analyzer.  A  linear  relationship  was  established  as  is 
shown  in  Fig.  BU.  The  graph,  Bit,  was  used  to  reduce  the  natural  fre¬ 
quencies  obtained  by  the  oscilloscope  to  the  corresponding  vibration 
analyzer  frequencies. 
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Bit  Calibration  of  Oscilloscope  versus  Vibration  Analyzer 
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